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Art Unit: 2121 

DETAILED ACTION 

1. The following is an initial Office Action upon examination of the above- 
identified application on the merits. Claims 1-11 are pending in this application. 

Priority 

2. Acknowledgment is made of applicant's claim for foreign priority based on an 
application filed in Canada on 3 September 2002. It is noted, however, that 
applicant has not filed a certified copy of the 2,400,580 application as required by 
35 U.S.C. 119(b). 

■ 

Information Oisclosure Statement 

3. The listing of ref evznczs in the specification is not a proper information 
disclosure statement. 37 CFR 1.98(b) requires a list of all patents, publications, or 
other information submitted for consideration by the Office, and MPEP 

§ 609.04(a) states, **the list may not be incorporated into the specification but 
must be submitted in a separate paper." Therefore, unless the references have 
been cited by the examiner on form PTO-892, they have not been considered. 
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Claim Objections 

4. Claim 4 is objected to because of the following informalities: 
system/method should be system. Appropriate correction is required. 

Claim Rejections - 35 USC § 112 

5. The following is a quotation of the second paragraph of 35 U.S.C. 112: 

The specification shall conclude with one or more claims particularly pointing out and distinctly 
claiming the subject matter which the applicant regards as his invention. 

6. Claim 1 is rejected under 35 U.S.C. 112, second paragraph, as being indefinite 
for failing to particularly point out and distinctly claim the subject matter which 
applicant regards as the invention. 

Regarding claim 1, the term "etc." vzx\6zvs the claim(s) indefinite because 
the claim(s) include(s) elements not actually disclosed (those encompassed by 
"etc."). thereby rendering the scope of the claim(s) unascertainable. See MPEP 
§ 2173.05(d). 

Allowable Subject Matter 

7. Claims 9-11 are allowed. 
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8. Claim 1 would be allowable If rewritten or amended to overcome the 
rejection(s) under 35 U.S.C. 112, 2nd paragraph, set forth in this Office action. 

9. Claim 2-8 would be allowable if rewritten to overcome the rejection(s) under 
35 U.S.C. 112, 2nd paragraph, set forth in this Office action and to include all of 
the limitations of the base claim and any intervening claims. 

10. The following is a statement of reasons for the indication of allowable 
subject matter: 

As per claim 1, the prior art of record taken alone or in combination fails to 
teach performing Loadf low computation at said nodes of the power system by a 
Super Super Decoupled computation in any of the Super Super Decoupled Loadf low 
methods or any of their hybrid combination or simple variants employing 
corresponding gain matrices derived from a super decoupled Jacobian matrix for 
real power with respect to angle and a super decoupled Jacobian matrix for 
reactive power with respect to voltage, and involving triangular factorization of 
said gain matrices and computing of discrepancy of real power and reactive power 
from specified values through such nodes, said computing including introducing 
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variables representing quotients of the transformed discrepancies from specified 
values of real and reactive power flowing in through each node divided by voltage, 
or square of the voltage in case of transformed real power mismatches in methods 
employing (l.theta., IV) iteration scheme, on each node, and using such variables to 
calculate values of angle and voltage for said transformed discrepancies from 
specified values of real and reactive power flowing in through each node, by using 
triangular factorization of said gain matrices for real and reactive power and 
restricting nodal transformation angle to maximum -48 degrees, applied to complex 
power injection in computing said transformed discrepancies from specified values 
of real and reactive power flowing in through each node. 

As per claim 9, the prior art of record taken alone or in combination fails to 
teach controlling the operation of the excitation element of at least one machine 
to produce or absorb the amount of reactive power indicated by any of the said 
Super Super decoupled models with respect to the set of specified parameters. 

Conc/us/of! 

11. The prior art made of record and not relied upon is considered pertinent to 
applicant's disclosure. 
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The following references are cited to further show the state of the art with 
respect to electrical power system regulation: 

USPN 7,096,175 B2 to Rehtanz et al. 
USPN 7.096,165 B2 to Pantenburg et al. 
USPN 6,754,597 B2 to Bertsch et al. 
USPN 6,690,175 B2 to Pinzon et al. 
USPN 6,313,752 Bl to Corrigan et al. 
USPN 5,610,834 to Schlueter 
USPN 5,566,085 to Marceau et al. 
USPN 5,305,174 to Morita et al. 
USPN 5,081,591 to Hanway et al. 
USPN 4,974.140 to Iba et al. 
USPN 4.868.410 to Nakamuro 
US Pub. No. 2007/0203658 Al to Patel 
WO 2004/023622 A2 to PATEL 
CA 2107388 A to PATEL 
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Art Unit: 2121 

Any inquiry concerning this communication or earlier communications from 
the examiner should be directed to Crystal J. Barnes Bullock whose telephone 
number is 571.272.3679. The examiner can normally be reached on Monday- Friday. 

If attempts to reach the examiner by telephone are unsuccessful, the 
examiners supervisor, Albert Decady can be reached on 571.272.3819. The fax 
phone number for the organization where this application or proceeding is assigned 
is 571-273-8300. 
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Information regarding the status of an application may be obtained from the 
Patent Application Information Retrieval (PAIR) system. Status information for 
published applications may be obtained from either Private PAIR or Public PAIR. 
Status information for unpublished applications is available through Private PAIR 
only. For more information about the PAIR system, see http://pair- 
direct.uspto.gov. Should you have questions on access to the Private PAIR system, 
contact the Electronic Business Center (EBC) at 866-217-9197 (toll-free). If you 
would like assistance from a USPTO Customer Service Representative or access to 
the automated information system, call 800-786-9199 (IN USA OR CANADA) or 
571-272-1000. 

/Crystal J. Barnes Bullock/ 
Primary Examiner, Art Unit 2121 
21 November 2008 
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(57) Abstract: Load-Flow computations are performed in 
reai-time operation/control and in on-line/off-line studies 
of electrical power systems. Three Load-Flow computation 
methods of the present invention are the best versions of many 
simple variants with almost similar performance. The Super 
Super Decoupled Loadflow (SSDL-YY) method and its many 
variants are characterized in limiting rotation angle applied to 
nodal real and reactive power mismatches to the maximum of 
-48 degrees instead of -36 degrees, replacing network shunt 
element -2b^Cos*p by [2(QSHpCos*p - PSHpSin*p)A^,2] or 
by [-b^Cos*p + (QSHpCos^p - PSHpSin$p)A^»2] and using die 
dividing term V* instead of V in the modified nodal real power 
residues [RP] in the prior art Fast Super Decoupled Loadflow 
(FSDL) method. The other two Super Super Decoupled 
Loadflow: BGX* version (SSDL-BGX') and X'GpvX' version 
(SSDL-X'GpvX*) are characterized in the use of simultaneous 
(IV, 19) iteration scheme thereby reducing the mismatch 
computation once compared to two mismatch computations 
in the prior art method employing successive (18,1V) iteration 
scheme. The invented methods are also characterized in the 
different definition of gain matrices as detailed in the steps of 
algorithm-2, aIgorithm-3 and algorithm-4 of carrying out of 
the inventions, (steps-cc, -dd, and ff in Rg.2; steps-ccc, -ddd, 
-fff, -ggg and -hhh in Fig.3; steps -dddd and -fiff in Rg. 4) 
leading to some speed-up of the invented methods. 
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SYSTEM OF SUPER SUPER DECOUPLED LO ADFLOW COMPUTATION FOR 

ELECTRICAL POWER SYSTEM 

TECHNICAL FIELD 

The preseat invention relates to methods of loadflow computation in power flow control 
and voltage control for an electrical power system. 

BACKGROUND ART AND MOTIVATION 

Utility/industrial power networks are composed of many power plants/generators 
interconnected through transmission/distribution lines to other loads and motors. Bach of 
these components of the power netwoik is protected against unhealthy (faulty, over/under 
voltage, over loaded) potentially damaging operating conditions. Such protection is 
automatic and operates without the consent of power network operator, and takes an 
unhealthy component out of service disconnecting it from the network. The time domain of 
operation of the protection is of the order of milliseconds. 

The purpose of a utility/industrial power netwoik is to meet the electricity demands of its 
various consumers 24-hours a day, 7-days a week while maintaining the quality of 
electricity supply. The quality of electricity supply means the consumer demands be met at 
specified (say + or - 5% tolerance) voltage and frequency levels without over loaded, 
under/over voltage operation of any of the power network components. The operation of a 
power network is different at different times due to changirig consumer demands and/or 
development of any faulty/contingency situation. Jn other words healthy opemting power 
network is constantiy subjected to small or large disturbances. These disturbances could be 
operator initiated, or initiated by security control functions and various optimization 
functions such as economic operation, minimization of losses etc., or caused by a 
fault/contingency incident 
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For example, a power network is operating healthy and meeting quality electricity needs of 
its consumers. A fault occurs on a liue or a transfoimer or a generator which faulty 
component gets isolated from the rest of the healthy network by virtue of the automatic 
operation of its protection. Such a disturbance would cause a change in the pattern of 
power flows in the network, which can cause over loading of one or more of the other 
components and/or over/under voltage at one or more nodes in the rest of the network This 
in turn can isolate one or more other components out of service by virtue of the operation 
of associated protection, which disturbance can trigger chain reaction disintegrating the 
power network 

Therefore, the most basic and integral part of all other functions (e.g. optimizations) in 
power network operation and control is security control. Security control means controlling 
power flows so that no component of the network is over loaded and controlling voltages 
such that ihc£Q is no over voltage or under voltage at any of the nodes in the network 
following a disturbance small or large. Security control functions (overload alleviation and 
over/under voltage alleviation) can be realized through one or combination of more 
controls ia the network. These involve control of power flow over tie liae coimecting other 
utility network, turbine steam/water input control to control real power generated by each 
generator, load shedding function curtails load demands of consnnfiers, excitation controls 
reactive power generated by individual generator which essentially controls generator 
temunal voltage, transformer taps control connected node voltage, switching in/out in 
capacitor/reactor banks controls reactive power at the connected node. Such overload and 
under/over voltage alleviation functions produce control amount changes in step-60 of 
Fig. 5. These control amount changes could be even optimized in case of simulation of all 
possible imaginable disturbances (outage of a line, loss of generation etc.) for corrective 
action stored and made readily available for acting upon in case the simulated disturbance 
actually occurs in the power network In fact simulation of all possible imaginable 
disturbances is tbe modem practice because corrective actions need be taken before the 
operation of individual protection of unhealthy component 
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Control of an electrical power system (Power-flow control, voltage control etc.) is 
performed according to tbie process flow diagram of fig.5. The various numbered steps in 
fig. 5 are explained in the following. 

Step-10: On-line readings of various real-time power flows, voltages, circuit breaker status 

(open/close) etc. are obtained 
Step-20: A control amount (Le. change in power injections, voltages etc.) is initially 

established and proposed 
Step-30: Various power flows, voltage magnitudes and angles, reactive power generations 

by generators and turns ratios of transformers in the power system are determined 

by performing loadflow computation, which incorporates established/proposed/set 

control adjustments 

Step-40: The results of Loadflow computation of step 30 are evaluated for any over loaded 

transmission lines and over/under voltages at different nodes in the power system 
Step-50: If the system state is good (no over loaded lines and no over/under voltages), the 

process branches to step 70, otherwise to 60 
Step-60: Changes the control amount initially set in step-20 or later set in the previous 

process cycle step-60 and returns to step-30 
Step-70: Actually implements the control amount correction to obtain 

securc/optimum/correct/acceptable operation of power system 

It is obvious that Loadflow computation is performed many times in real-time operation 
and control environment and, therefore, high-speed (effident) Loadflow computation is 
necessary to provide corrective control in the changing power system conditions including 
an outage or failure. Moreovo:, the loadflow compntations must be highly reliable to 
yield converged solution under wide range of system operating conditions and 
network parameters. Failure to yield converged loaflow solution creates blind spot as to 
what exacfly could be h^pening in the network leading to potentially damaging 
operational and ccmtrol decisions/actions in coital-intensive power utihties. 
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The embodiment of the present invention, the most efficient and reliable loadflow 
computations, as described in the above steps and in Fig. 5 is very general and elaborate. 
The control of voltage magnitude within reactive power generation c^abilities of electrical 
machines (generators, synchronous motors, capadtor/inductor banks) and within operating 
ranges of transfonner taps is normally integral part of Loadflow computation as described 
in '*LTC Transformers and MVAR violations in the Fast Decoupled Loadflow, IEEE PAS- 
101, No.9, PP. 3328-3332." If under/over voltage still exists in the results of Loadflow 
computations, other control actions are taken in step-60 in the above and in Fig.5. For 
example, xmder voltage can be alleviated by shedding some of the load connected 

However, the simplest embodiment of the efficient and reliable system and method of 
loadflow computations is where only voltage magnitudes are controlled by controlling 
reactive power generation of generators and motors, switching in/out in capacitor/inductor 
banks and transformer taps. Of course, such control is possible only within reactive power 
capabilities of machines and capacitor/reactor banks, and within operating ranges of 
transformer taps. This is the case of a network in which the.real power assignments have 
already been fixed and in which steps-50 and -60 in the above and in Fig.5 are absent Once 
loadflow computations are finished, the Loadflow solution includes indications of reactive 
power generation at generator nodes and at the nodes of the capadtor/inductor banks, and 
indications of transfctmer tap settings. Based on the known reactive power capability 
characteristics of the individual machines, the determined reactive power values are used to 
adjust the excitation current to each machine, or at least each machine presently under 
reactive power, or VAR, control, to establish the desired reactive power set points. The 
transformer taps are set in accordance with the tap setting indications produced by the 
Loadflow computation system. 

This procedure can be employed either on-line or off-line. In off-Une operation, the user 
can simulate and experiment with various sets of operating conditions and determine 
reactive power generation and transformer tap settings requirements. A general-purpose 
computer can implement the entire system. For on-line operation, the loadflow computation 
system is provided with data identifying the current real and reactive power assigmnents 

4 
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and transformer transformation ratios, the present status of all switches and circuit breakers 
in the network and machine characteristic curves in steps-10 and -20 in the above and in 
Fig. 5, and blocks 10, 12, 14, 20, 30, 40, 42, 44, 50, 52 in Fig 6. Based on this information, 
a model of the system based on gain matrices of any of the invented or prior art Loadflow 
computation methods provide the values for the corresponding node voltages, reactive 
power set points for each machine and the transformation ratio and tap changer position for 
each transformer. 

m 

The present invention relates to control of utility/industrial powCT networks of the types 
including plurality of power plants/generators and one or more motors/loads, and 
connected to other external utility. In the utility/industrial systems of this type it is the usual 
practice to adjust the real and reactive power produced by each generator and each of the 
other sources (synchronous condensers, capacitor/inductor banks) in order to optimize the 
real and reactive power generation assignments of the system. Healthy (secure) operation 
of the network can be shifted to optimized operation through corrective control 
(disturbance) produced by optimization functions without violation of security constraints. 
This is referred to as security constrained optimization of operatioiL Such an optimization 
is described in the United States Patent Number. 5,081,591 dated Jan. 13, 1992 
(Optimizmg Reactive Power Distribution in an Industrial Power Networic) where the 
present invmtion can be embodied by replacing the block nos. 56 and 66 by a block of 
constant matrices [YG] and [YV], and replacing blocks of 'TExercise Newton-Raphson 
Algorithm" by blocks of 'Exercise Fast Super Decoupled Algorithm" or '^Execise Super 
Super Decoupled Algorithm" in place of blocks 58 and 68. 

DISCLOSURE OF THE INVENTION 

This invention relates to steady-state power network computation referred to as Loadflow 
or Power-Flow. Loadflow computations are performed as a step in real-time 
operation/control and in on-line/off-line studies of Electrical Power Systems. The present 
invention involves three-methods. These 3-methods are the best versions of many simple 
variants with almost similar perfomciance. Simple variants include any possible hybrid 

5 
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combination of these 3 -methods and unsymmetiical definitions of [Y6] in SSDL-method. 
Among these 3-methods, their variants and all other known methods. Super Super 
Decoupled Loadflow (SSDL-YY) is the simplest, easiest to implement and overall best 
in performance (reliability of convergence and efficiency of computations). 

BRIEF DESCRIPTION OF DRAWINGS 

Fig. 1 is a flow-chart of the prior art Fast Super Decoupled Loadflow computation method 
Fig. 2 is a flow-chart embodiment of the invented Super Super Decoupled Loadflow 

computation method of version SSDL-YY 
Fig. 3 is a flow-chart embodiment of the invented Super Super Decoupled Loadflow: 

SSDL-BGX\ -BGY, and -BGX versions 
Fig. 4 is a flow-chart embodiment of the mvented Super Super Decoupled Loadflow: 

SSDL-X'GpvX', SSDL-YGpyY, and SSDL- XGpvX versions 
Fig. 5 is a flow-chart of the overall controlling method for an electrical power system 

iavolving Loadflow computation as a step which can be executed using one of the 

Loadflow computation methods of Figs. 1, 2, 3, 4, other variations described, their 

hybrid combination and/or their simple variants 
Fig. 6 is a flow-chart of the simple special case of voltage control in overall controlling 

method of Fig. S for an electrical power system 

Symbols 

* 

The prior art and inventions will now be described using the following symbols: 

Ypq = Gpq + jBpq : (p-q) th element of nodal admittance matrix without shunts 

y = gp + jbp : total shunt admittance at any node-p 

Vp = ^ + jfp = VpZ0p : complex voltage of any node-p 

Vs = + jfc = VsZGs : complex slack-node voltage 

A0p, AVp : voltage angle, magnitude corrections 

Acp, : real, imaginary components of voltage corrections 
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I*p jQp ■ i^odal injected power, calculated 

APp + j AQp : nodal power residue (mismatch) 

RPp + jRQp : modified nodal power residue 

PSH^ + jQSHp : net nodal injected power, scheduled 

Op : rotation angle 

m : number of PQ-nodes 

k : number of PV-nodes 

n=m+k+l : total number of nodes 

q>p : q is the node adjacent to node-p excluding the case of q=p 

[ ] : indicates enclosed variable symbol to be a vector or a matrix 

LRA : Limiting Rotation Angle 



PQ-node : load-node (Real-Power-P and Reactive-Power-Q are specified) 
PV-node : generator-node (Real-Power-P and Voltage-Magnitude-V are specified) 

Decoupled Loadflow 

A class of decoupled Loadflow methods involves a system of equations for the separate 
calculation of voltage angle and voltage magnitude corrections. Each decoupled method 
conq)rises a system of equations (1) and (2) differing in the definition of elements of [RP], 
[RQ], and [YG] and [YV]. 

4 

[RP] = [Ye] [Ae] (1) 

[RQ] = [YV] [AV] (2) 
Successive (19, IV) Iteration Scheme 

In this scheme (1) and (2) are solved alternately with intermediate updating. Each iteration 
involves one calculation of [RP] and [A6] to update [6] and then one calculation of [RQ] 
and [AV] to update [V]. The sequence of relations (3) to (6) depicts the scheme. 
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[A8] = [Ye] * [RP] 



(3) 



[6] = [6] + [Ae] 



(4) 



[AV] = [YV] [RQ] 



(5) 



tV] = [V] + [AV] 



(6) 



The scheme involves solution of system of equations (1) and (2) in an iterative manner 
depicted in the sequence of relations (3) to (6). This scheme requires mismatch calculation 
for each half iteration; because \BJ?] and [RQ] are calculated always using the most recent 
voltage values and it is block Gauss-Seidal approach. The scheme is block successive, 
which imparts ino-eased stabihty to the solution process. This in turn improves 
convergence and increases the reliability of obtaining solution. 

PRIOR ART: FAST SUPER DECOUPLED LOADFLOW METHO 
(Referen.ces-3, 6, 7) 

Fast Snper Decoupled Loadflow (FSDL) Method 

RPp = (APpCosOp + AQ,SinOp)A^p -for PQ-nodes (7) 
RQ, = (AQpCosOp - APpSinOp )A^p -for PQ-nodes (8) 
CosOp = Absolute (Bpp / /(Gpp* + Bpp^)) ^ Cos (-36°) (9) 




(10) 



RPp = APp/(KpVp) 



-for PV-nodes 



(11) 
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YGpq- 



-Ypq -for branch r/x ratio ^ 2.0 

- (Bpq + 0.9(Ypq-Bpq)) -foT branch r/x ratio > 2.0 

^pq -for branches connected between two PV-nodes or 

a PV-node and the slack-node (12) 



YVpq = \-Ypq -for branch r/x ratio ^ 2.0 

L (Bpq + 0.9(Ypq-Bpq)) -foT branch r/x ratio > 2.0 (13) 



Yepp = S-Yepq and YVpp = -2bp' + S-YVpq (14) 

q>p q>p 



bp' = bpCosOp or bp* = bp (15) 

« 

= Absolute (Bpp/Y0pp) ( 1 6) 

t 

% 

Branch admittance magnitude ia (12) and (13) is of the same algebraic sign as its 
susceptance. Elements of the two gain matrices differ in that diagonal elements of [YV] 
additionally contain the b' values given by relation (IS) and in respect of elements 
corresponding to branches connected between two PV-nodes or a PV-node and the slack- 
node. Relations (9) and (10) with inequality sign implies that nodal rotation angles are 
restricted to maximum of -36 degrees. The method consists of relations (3) to (16). In two 
simple variations of the FSDL method, one is to make YVpq=Y6pq and the other is to make 
YGpq-YVpq. Kp is restricted to the minimum value of 0.75 detemiined experimentaUy, and 
it is system independent However it can be tuned for the best possible convergence for any 
given system. 

This prior art method involves solution of system of equations (1) and (2) in an iterative 
manner depicted in sequence of relations (3) to (6). Prior art method is embodied in 
algorithm-1, and in the flow-chart of fig. 1. 



Computation steps of FSDL method (AIgorithm-1): 
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a. Read system data and assign an initial approxiniate solution. If better solution 
estimate is not available, set voltage magnitude and angle of all nodes equal to those 
of the slack-node. This is referred to as the slack-start 

b. Fonn nodal admittance matrix, and Initialize iteration count ITRP = ITRQ= r = 0 

a Compute Cosine and Sine of nodal rotation angles using relations (9) and (10), and 
store them. If they, respectively, are less than the Cosine and Sine of -36 degrees, 
equate them, respectively, to those of -36 degrees. 

d. Form (m+k) x (m+k) size matrices [Y9] and [YVJ of (1) and (2) respectively each 
in a compact storage exploiting sparsity. The matrices are formed using relations 
(12), (13), (14), and (15). In [YV] matrix, replace diagonal elements corresponding 
to PV-nodes by very large value (say, 10.0**10). In' case [YV] is of dimension (m x 
m), this is not required to be performed. Factorize [Y9] and [YV] using the same 
ordering of nodes regardless of node-types and store them using the same 
indexing and addressing information. In case [YV] is of dimension (m x m), it is 
factorized using different ordering than that of [Y9]. 

e. Compute residues [AP] (PQ- and PV-nodes) and [AQ] (at only PQ-nodes). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step. 

f. Compute the vector of modified residues [RP] using (7) for PQ-nodes, and using 
(1 1) and (16) for PV-nodes. 

g. Solve (3) for [A9] and update voltage angles using, [9] = [9] + [A9]. 

h. Set voltage magnitudes of PV-nodes equal to the specified values, and Increment 
the iteration count ITRP=ITRP+1 and r=(lTRP+ITRQ)/2, 

L Compute residues [AP] (PQ- and PV-nodes) and [AQ] (at PQ-nodes only). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step. 

j. Compute the vector of modified residues [RQ] using (8) for only PQ-nodes. 

k. Solve (5) for [AV] and update PQ-node magnitudes using [V] = [V] + [AV]. While 
solving equation (5), skip all the rows and columns corresponding to PV-nodes. 
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1. Increment the iteration count ITRQ=ITRQ+1 and i=(ITRP+ITRQ)/2, and Proceed 
to step (e) 

m. Calculate line flows and output the desired results 



BSrVENTED SUPER SUPER DECOUPLED LOADFLOW METHODS 



Super Super Decoupled LoadQow: X'X'-version (SSDL-X'X') 

The general method, in successive (16, IV) iteration scheme represented by sequence of 
relations (3) to (6), can be realized as SSDL-X'X', from which manifested are many 
versions. The elements of [RP], [RQ], [YG] aad [YV] are defined by (17) to (29). 

RPp = [APp' + (GppV Bpp') AQp'] / Vp^ -for PQ-nodes (17) 

RQp = [AQp'-(Gpp'/Bpp')APp']/Vp -for PQ-nodes (18) 

RPp= [APp / (Kp*Vp^)] -forPV-nodes (19) 

Y%^= -l/Xpq' aad YVpq=-l/Xpq' (20) 



Yepp = S-Yep, and YVpp = bp' + 2-YVpq (21) 



Whore, 

bp' = -2bpCos<E>p or 
bp' = -b„CosOp + [QSI^' - ((^ VBpp') PSHp'] / V,^ or 

bp' = 2[QSHp' - (GppV Bpp') PSHj,'] / V,^ (22) 
APp' = APpCosOp + AQpSin^p -for PQ-nodes (23) 

AQp' = AQpCosOp - APpSm<^ -for PQ-nodes (24) 

11 



SUBSTITUTE SHEET (RULE 26) 



wo 2004/023622 PCT/CA2003/001312 

PSHp' = PSHpCosOp + QSHpSinOp -for PQ-nodes (25) 

QSHp' = QSHpCosOp - PSHpSin^^p -for PQ-nodes (26) 

CosOp = Absolute [Bpp / v (Gpp^ + Bpp^)] ^ Cos (any angle from 0 to -90 degrees) (27) 

SinOp = -Absolute [Gpp / V (Gpp^ + Bpp^)] ^ Sin (any angle from 0 to -90 degrees) (28) 

Kp = Absolute (BppA^Opp) (29) 



The factor of (29) is initially restricted to the TniniTnum of 0.75 detemiined 
e2q)erimentally; however its restriction is lowered to the minimum value of 0.6 when its 
average over all PV-nodes is less than 0.6. This factor is system and method independent. 
However it can be tuned for the best possible convergence for any given system. This 
statement is valid when the factor is applied in the manner of equation (19) in all the 
methods derived in the following from the most general method SSDL-X'X\ 

The definition of YBpq in (20) is simplified because it does not explicitly state that it always 
takes the value of -Bpq for a branch connected between two PV-nodes or a PV-node and the 
slack-node. This fact should be understood implied in all the definitions of Y6pq in this 
document 

However, a whole new class of methods, corresponding to all those derived in the 
following and prior art, results when the factor Kp is used as a multiplier in the definition of 
RPp at PQ-nodes as in (30) instead of divider in RPp at PV-nodes as given in (19). This 
will cause changes only in relations (17), (19), and (20) as given in (30), (31), and (32). 

RPp = {[APp' + (Gpp' /Bpp') AQp'] / Vp2 }* Kp -for PQ-nodes (30) 

KPp= APp/Vp2 -for PV-nodes (31) 
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Y0pq= -Bpq and YVp<i = (32) 



The best performance of methods of this new class has been realized when tihie factor E^, 
appUed in a manner of relation (30) leading to changes as in (30) to (32), is unrestricted. 
That means it can take any value as given by relation (29). 

Super Super Decoupled Loadflow: YY-version (SSDL-YY) 



If unrestricted rotation is applied and transformed susceptance is taken as admittance value 
with the same algebraic sign and transformed conductance is assumed zero (reference-6), 
the SSDLrX'X' method reduces to SSDL-YY. Though, this method is not very sensitive to 
the restriction applied to nodal rotation angles, SSDL-YY presented here is the best 
possible experimentally arrived at method However, it gives closely similar performance 
over wide range of restriction applied to the nodal rotation angles (say from -36 to -90 
degrees). 

RPp = APp7Vp^ and RQp==AQpVVp -forPQ-nodes (33) 

RPp = APp / (KpVp^) -for PV-nodes (34) 

Y6pq = \ -Yp(i -for branch r/x ratio < 3.0 

- (Bpq + 0.9(Ypq-Bpq)) -foT branch r/x ratio > 3.0 

^-Bpq -for branches connected between two PV-nodes or 

a PV-node and the slack-node (35) 



YVpq = 1 -Ypq -for branch r/x ratio <. 3.0 

- (Bpq + 0.9(Ypq-Bpq)) -for branch r/x ratio > 3.0 (36) 

Yepp^Z-YBpq and YVpp = bp* + S-YVpq (37) 

> 

bp' = (QSHpVVs^)-bpCosOp or bp' = 2QSHp V V,» (38) 
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where, APp', AQp' , QSHp' , and are defined in relations (23) to (29). However, nodal 
rotation angles in relations (27) and (28) be restricted to the maximum of -48 degrees for 
this method, determined experimentally for the best possible convergence firom non 
linearity considerations. 

In case of systems of only PQ-nodes and without any PV-nodes, equations (35) and 
(36) simply be taken as YOpq = -Ypq and YVpq = -Ypq • The factor Kp of (29) is initiany 
restricted to the minliniini of 0.75 determined experimentally; however its restriction 
is lowered to the mip™™ value of 0.6 when its average over all PV nodes is less than 
0.6. This factor is system independent HowevCT it can be tuned for the best possible 
convergmce for any given systeoL 

Branch admittance magnitude in (35) and (36) is of the same algebraic sign as its 
susceptance. Elements of the two gain matrices differ in that diagonal elements of [YV] 
additionally contain the b' values given by relations (37) and (38) and in respect of 
elements corresponding to branches connected between two PV-nodes or a PV-node and 
the slack-node. Relations (27) and (28) with inequality sign implies that nodal rotation 
angles are restricted to maximum of -48 degrees for SSDL-YY. The method consists of 
relation's (3) to (6), (33) to (38), and (23) to (29). In two simple variations of the SSDL- 
YY metiiod, one is to make YVpq=Yepq and the other is to make Yepq=YVpq. 

SSDL-YY method implements successive (10, IV) iteration scheme and is embodied in 
algorithm-2, and in flow-chart of fig.2 where double lettered steps are characteristic steps 
of the SSDL-YY method and are different than those of the prior art FSDL method 

Computation steps of SSDL-YY method (Algorifhm-2): 

a. Read system data and assign an initial approximate solution. If better solution 
estimate is not available, set voltage magnitude and angle of all nodes equal to those 
of the slack-node. This is referred to as the slack-start 

b. Form nodal adnaittance matrix, and Initialize iteration count ITRP = ITRQ= r = 0 
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cc Compute Cosine and Sine of nodal rotation angles using relations (27) and (28), and 
store thenL If ihey, respectively, are less than the Cosine and Sine of -48 degrees, 
equate them, respectively, to those of -48 degrees. 

dd. Form (m+k) x (m+k) size matrices [Y8] and [YV] of (1) and (2) respectively each 
in a compact storage exploiting sparsity. The matrices are formed using relations 
(35), (36), (37), and (38). In [YV] matrix, replace diagonal elements corresponding 
to PV-nodes by very large value (say, 10,0**10). In case [YV] is of dimension (m x 
m), this is not required to be performed. Factorize [Y9] and [YV] using the same 
ordering of nodes regardless of node-types and store them using the same 
indexing and addressing information* in case [YV] is of dimension (m x m), it is 
factorized using different ordering than that of [Y9]. 

e. Compute residues AP (PQ- and PV-nodes) and AQ (at only PQ-nodes). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step. 

ff. Compute the vector of modified residues [RP] as in (33) for PQ-nodes, and using 
(34) and (29) forPV-nodes. 

g. Solve (3) for [AO] and update voltage angles using, [6] = [6] + [AO]. 

h. Set voltage magnitudes of PV-nodes equal to the specified values, aud Increment 
the iteration count ITRP=ITRP+1 and i=(ITRP+ITRQ)/2. 

i. Compute residues [AP] (PQ- and PV-nodes) and [AQ] (at PQ-nodes only). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step. 

j . Compute the vector of modified residues [RQ] as in (33) for only PQ-nodes. 

k Solve (5) for [AV] and update PQ-node magnitudes using [V] = [V] + [AV], While 

solving equation (5), skip all the rows and coluimis corresponding to PV-nodes. 
1. Increment the iteration count ITRQ=rrRQf 1 and i=(ITRP+ITRQ)/2, and Proceed 

to step (e) 

m Calculate line flows and oulput the desired results 



Super Super Decoupled Loadflow: XX-version (SSDL-XX) 

If no (zero) rotation is applied, the SSDLrX'X' method reduces to. SSDL-XX, which is the 
simplest form of SSDL-X'X\ The SSDL-XX method comprises relations (3) to (6), (39) to 
(45), and (29). 
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RPp = [APp + (Gpp / Bpp) AQp] / Vp^ -for PQ-nodes (39) 

RQp = [AQp-(Gpp/Bpp)APp]/Vp -for PQ-nodes (40) 

RPp = APp / (KpVp^) -for PV-nodes (41) 



YGpq- 



- 1 .0/Xpq -for all other branches 



^-Bpq -for branches connected between two PV-nodes or 

a PV-node and the slack-node (42) 



YVpq= -1.0/Xpq -for aU branches (43) 



Yepp = E-Yepq and YVpp = bp» + 2-YVpq (44) 

<f>p q>p 



bp* = -2bp or 
V = -bp + [QSHp - (Gpp / Bpp) PSHp]A^s^ or 

bp* = 2[QSHp - (Gpp / Bpp) PSHp]A^s^ (45) 

where, Kp is defined in relation (29). This is the simplest method with very good 
perfonnance for distribution networks in absence of PV-nodes (for systems containing only 
PQ-nodes). The large value of the difference [(1/X)-B], particularly for high R/X ratios 
branches connected to PV-nodes, creates modeling error when PV-nodes are present in a 
systenx 



Super Super Decoupled Loadflow: BX-version (SSDLrBX) 



If super decoupling is applied only to QV-sub problem, the SSDL-XX method reduces to 
SSDL-BX, which makes it perform better for systems containing PV-nodes. The SSDL- 
BX method comprises relations (3) to (6), (46) to (48), (44) and (45). This method can be 
referred to as Advanced BX-Fast Decoupled Loadflow. 
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(46) 



RQp = [AQp - (Gpp / Bpp) APp] / Vp -for PQ-nodes (47) 

Yep,= -Bpq and YVp, = -l/Xpq (48) 

It should be noted that Amerongen's General-purpose Fast Decoupled Loadflow method of 
reference-5 has turned out to be an approximation of this method. The approximation 
involved is only in relation (47), However, numerical performance is found to be only 
slightly better but more reliable than that of the Amerongen's method 

Super Super Decoupled Loadflow: X'B'-versiou (SSDL-X'B') 

RPp = [APp' + (Gpp V Bpp') AQp*] / Vp^ -for PQ-nodes (49) 

RQp= AQpVVp -for PQ-nodes (50) 

RPp= [APp/(I^*Vp^)] -for PV-nodes (51) 

Yepq=-1/Xp,' and YVp, = -Bp,' (52) 

Yepp = Z-Yepq and YVpp = V + I^-YVpq (53) 

Where, bp' = -2bpCos$p or 

bp' = -bpCosOp + QSHpV V/ or 

bp' = 2QSHpV V,^ (54) 

Where, APp', AQp', PSHp', QSHp', CosOp, SinOp, Kp are defined in (23)10(29).^ method 
consists of relations (3) to (6), (49) to (54), and (23) to (29). Best performance of fliis 
method could be achieved by restricting $p in (27) and (28) to less than or equal to -48°. 
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Super Super Decoupled Loadflow: YB'-version (SSDLrYB') 



The relation (49) in SSDL-X'B' implies unrestricted Op is applied and it can take values up 
to -90 degrees. Therefore, (49) can be modified to (55) with consequent modification of 
(52) into (56). 

RPp = [APp*Absolute [Bpp / v (Gpp^ + Bpp^)] + AQp* [-Absolute [Bpp / v (Gpp^ + Bpp^)]] / Vp^ 

-for PQ-nodes (55) 

Y0pq= -Ypq and YVpq = -Bpq' (56) 

This method consists of relations (3) to (6), (55), (50), (51), (56), (53) and (54), and (23) to 
(29). Best performance of this method could be achieved by restricting Op in (27) and (28) 
to less than or equal to -48 degrees. Where, APp' , AQp' , PSHp' , QSHp' , CosOp , SinOp , 



are defined in (23) to (29). 
Super Super Decoupled Loadflow: B'X'-versiou (SSDL-B'X') 

RPp = APpV Vp^ -for PQ-nodes (57) 

RQp = [AQp' - (Gpp'/ BppO APp'] / Vp -for PQ-nodes (58) 

RPp = [APp / (Kp*Vp^)] -for PV-nodes (59) 

YSp, = -Bp,' and ' YVp, = -l/X,,' (60) 

Yepp = Z-Yep, aad YVpp = bp' + S-YVpq (61) 

Where, bp' = -2bi,CosOp or 

V = -bpCosOp + [QSHp' - (Gpp VBpp') PSHp'yVs^ or 
bp' = 2[QSHp' - (GppVBpp') PSHp'yV,^ (62) 
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Where, APp', AQp', PSHp', QSHp', CosOp, SinOp, Kp are defined in (23)to (29).This 
method consists of relations (3) to (6), (57) to (62), and (23) to (29). Best perfoimance of 
this method could be achieved by restricting Op in (27) and (28) to less than equal to -48°. 

Super Super Decoupled Loadflow: B' Y-version (SSDL- B'Y) 

The realtion (58) m SSDL-B'X' implies unrestricted Op is applied and it can take values up 
to -90 degrees. Therefore, (58) can be modified to (63) with consequent modification of 
(60) into (64). 

RQp = [AQp*Absolute [Bpp / v (Gpp^ + Bpp^)] - APp* [-Absolute [Bpp / v (Gpp^ + Bpp^)]] / Vp^ 

-for PQ-nodes (63) 
Y9pq=-Bpq' and YVpq = -Ypq (64) 

This method consists of relations (3) to (6), (57), (63), (59), (64), (61) and (62), and (23) to 
(29). Best performance of this method could be achieved by restricting Op in (27) and (28) 
to less than or equal to -48 degrees. Where, APp' , AQp' , PSHp' , QSHp' , CosOp , SinOp , 
are defined in (23) to (29). 



Simultaneous (lY, 16) Iteration Scheme 

An ideal to be approached for the decoupled Loadflow methods is the constant matrix 
Loadflow of refermce-6 referred in this document as BGGB-mefhod. In an attempt to 
imitate it, a decoupled class of methods with simultaneous (IV, 16) iteration scheme 
depicted by sequence of relations (65) to (69) is invented This scheme involves only one 
mismatch calculation in an iteration. The correction vector is calculated in two separate 
parts without intermediate updating. Each iteration involves one calculation of [RQ], [AY], 
and [RP], [A8] to update [V] and [6]. 

[AY] = [YY] [RQ] (65) 
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[RP] = [APATI - [G] [AV] (66) 

[AG] = [Ye] [RP] (67) 

t 

[6] = [0] + [A9] (68) 

[V] = [V] + [AV] (69) 



In this invented class, each method differs only in the delBnition of elements of [RQ] and 
[YV]. The accuracy of methods depends only on the accuracy of calculation of [AV]. The 
greater the angular spread of branches terminating at PQ-nodes, the greater the inaccuracy 
in the calculation of [AV]. 

Super Super Decoupled Loadflow: BGX'-version (SSDL-BGX*) 

Numerical perfomiance could further be improved by cnganizing the solution in a 
simultaneous (IV, 19) iteration schone represented by sequence of relations (65) to (69). 
The elements of [RP], [RQ], [YB] and [YV] are defined by (70) to (74). 

(70) 

(71) 
(72) 

Yepp = E-Yepq and YVpp = bp' + E- YVp, (73) 

bp' = -IbpCosOp or 
bp' = - bpCosOp + [QSI^' - (Gpp VBpp') PSHp']A^,^ or 

bp' = 2[QSHp' - (GppVBppO PSHp']/ V,^ (74) 
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Where, APp' , AQp' , PSHp' , QSHp' , CosOp , SinOp are defined in (23) to (28).The SSDL- 
BGX' method comprises relations (65) to (74), and (23) to (28). Best possible convergence 
could be achieved by restricting rotations (Op) in the range (-10° to -20°) in relations (27) 
and (28). The method is embodied in algorithm-3 and in the flow-chart of Fig.3. 

Super Super Decoupled Loadflow: BGY-version (SSDL-BGY) 

If unrestricted rotation is applied and transformed susceptance is taken as admittance 
values and transformed conductance is assumed zero (reference-6), the SSDI^BGX' 
method reduces to SSDL-BGY as defined by relations (75) to (79). 

RQp = AQp Wp = (AQpCosOp - APpSinOp)A^p -for PQ-nodes (75) 

m 

RPp = (APpA^p) - SGpqAVq -for all nodes (76) 

Y0pq=-Bpq and YVpq = -Ypq (77) 

(78) 
(79) 

The SSDL-BGY method comprises relations (65) to (69), and (75) to (79). It is the special 
case of the SSDL-BGX' method 

Super Super Decoupled Loadflow: BGX-versiou (SSDL-BGX) 

If no (zero) rotation is appHed, the SSDL-BGX' method reduces to SSDL-BGX as defimed 
by relations (80) to (84). 
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RQp = [ AQp - (Gpp / Bpp)APp] / Vp -for PQ-nodes (80) 



m 



RPp= (APpA^p)-SGp,AVq -for aU nodes (8i) 



Yepq = -Bp, and YVpq = -l/Xp<, (82) 

Y6pp = S-Yep, and YVpp = V + S-YVp, (83) 

<PV q>p ^ ^ 

V = "2bpCosOp or 

V = -bpCos^p + [QSHp ~ (Gpp/Bpp)PSHp]A^s^ or 

V = 2[QSHp- (Gpp/Bpp)PSHp]A^s' (84) 

The SSDL-BGX method comprises relations (65) to (69), and (80) to (84). It is again the 
special case of the SSDL-BGX' method 



Compntation steps of SSDL-BGX% SSDL-BGY and SSDLrBGX methods 
(AIgorithm-3): 



a. Read system data and assign an initial approximate solution. If better solution 
estimate is not available, set voltage magnitude and angle of all nodes equal to those 
of the slack-node. This is referred to as the slack-start 

b. Form nodal admittance matrix, and Initialize iteration count ITR = 0. 

ccc. Compute Sine and Cosine of nodal rotation angles using relations (28) and (27), and 
store them. If they, respectively, are less than the Sine and Cosine of any angle set 
(say in the range -10 to -20 degrees), equate them, respectively, to those of the same 
angle in the range -10 to -20 degrees. In case of zero rotation. Sine and Cosine 
value vectors are not required to be stored 

ddd. Form (m+k) x (m+k) size matrices [YO] and [YV] of (1) and (2) respectively each 
in a compact storage exploiting sparsity 

1) In case of SSDL-BGX'-method, the matrices are formed using relations (72), 
(73), and (74) 
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2) In case of SSDL-BGY-method, the matrices are formed lising relations (77), 
(78), and (79) 

3) In case of SSDI/-BGX-niethod, the matrices are fonned using reladons (82), 
(83), and (84) 

In [YV] matrix, replace diagonal elements conesponding to PV-nodes by very large 
value (say, 10.0**10). In case [YV] is of dimension (m x m), this is not required to 
be performed. Factorize [Y6] and [YV] using the same ordering of nodes regardless 
of node-types and store them using the same indexing and addressing 
information. In case [YV] is of dimension (m x m), it is factorized using different 
ordering than that of [Y0]. 

e. Compute residues AP (PQ- and PV-nodes) and AQ (at only PQ-nodes). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step, 

ftt Compute the vector of modified residues [RQ] using (70) in case of SSDL-BGX' , 
using (75) in case of SSDL-BGY, and using (80) in case of SSDL-BGX method for 
only PQ-nodes, Solve (65) for [AVj. While solving equation (65), skip all the rows 
and columns corresponding to PV-nodes. Compute the vector of modified residues 
[RP] using (71) or (76) or (81). Solve (67) for [A9]. 

ggg- Update voltage angles using, [6] = [8] + [AO], and update PQ-node voltage 
magnitudes using [V] = [V] + [AVJ. 

hhh. Set voltage magnitudes of PV-nodes equal to the specified values, and Incremwt 
the iteration count ITR=ITR+1 , and proceed to step (e). 

m. Calculate line flows and output the desired results 

Triple lettered steps are characteristic steps of algorithm-3. The SSDL-BGX*, SSDL-BGY 
and SSDL-BGX methods differ only in steps-ccc and -ddd defining gain matrices, and 
step-fiff for calculating [RP] and [RQ]. Fig.3 is the flow-chart embodiment of algorithm-3. 

* < 

Super Super Decoupled Loadflow: X'GpvX'-version (SSDL- X'GpvX') 

Nimierical performance could also be improved by oiganizing the solution in a 
simultaneous (IV, 16) iteration schane represented by sequence of relatioiis (65) to (69). 
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The elements of [RP], [RQ], [Ye] and [YV] for this method are defined by (85) to (91). 



RPp = {[APp' + ((^'/Bpp')AQp'] / Vp^} - (gp'AVp) -for PQ-nodes (85) 



RQp = [AQp'-(Gpp/Bpp')APp']/Vp -for PQ-nodes (86) 

RPp= [(APpA^p^)-£GpqAVJ/Kp -forPV-nodes (87) 

Y8p, = -1/Xpq' and YVp, = -l/Xp,' (88) 



Y0pp = Z-Yepq and YVpp = V + S-YVp, (89) 



bp' = -2l^CosOp or 
bp' = -bpCosOp + [QSHp' - (Gpp V Bpp')PSHp']/Vg^ or 

bp' = 2[QSHp' - (Gpp VBpp')PSHp'yVs^ (90) 



gp' = 2bpSind>p or 
gp* = bpSinOp + [PSI^,' + (Gpp 7 BppOQSB^'yV.^ or 

gp' = 2[PSHp' + (Gpp'/Bpp')QSHp'yVs^ (91) 

Where, APp' , AQp' , PSHp' , QSE^' , CosOp , SinOp , are defined in (23) to (29) Again, 
if unrestricted rotation is applied and fransfonned susceptance is taken as admittance values 
and transformed conductance is assumed zero (reference-6), the SSDL-X'G^X' method 
reduces to SSDL-YGpvY. if no (zero) rotation is applied, the SSDL-X'GpvX' method 
reduces to SSDL-XGpvX. The SSDL-X'GpvX' method comprises relations (65) to (69), 
(85) to (91), and (23) to (29). It is embodied in algorithm-4 and in the flow-chart of Fig.4. 

Computation steps of SSDL-X'Gp,X', SSDIr-YGpvY and SSDL-XGt„rX methods 
(Algorithm-4): 
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a Read system data and assign an initial approxirDate solution. If better solution 
estbtnate is not available, set voltage magnitude and angle of all nodes equal to those 
of the slack-node. This is referred to as the slack-start 

b. Form nodal admittance matrix, and Initialize iteration count rrR= 0. 

ccc Compute Sine and Cosine of nodal rotation angles using relations (28) and (27), 
store them. If they, respectively, are less than the Sine and Cosine of any angle set 
(say 0 to -90 degrees), equate them, respectively, to those of the same angle in the 
range 0 to -90 degrees. In case of zero rotation, Sine and Cosine vectors are not 
required to be stored. 

dddd. Form (m+k) x (m+k) size matrices [YO] and [YV] of (1) and (2) respectively each 
in a compact storage exploiting sparsity using relations (88), (89), and (90). 
In [YV] matrix, replace diagonal elements corresponding to PV-nodes by very large 
value (say, 10.0**10). lii case [YV] is of dimension (m x m), this is not required to 
be performed. Factorize [YG] and [YV] using the same ordering of nodes regardless 
of node-types and store them using the same indexing and addressing 
information. In case [YV] is of dimension (m x m), it is factorized using different 
ordering than that of [Y9]. 

e. Compute residues AP (PQ- and PV-nodes) and AQ (at only PQ-nodes). If all are 
less than the tolerance (e), proceed to step (m). Otherwise follow the next step. 

mL Compute [RQ] usmg (86) for only PQ-nodes. Solve (65) for [AV]. While solving 
equation (65), skip all the rows and columns corresponding to PV-nodes. Compute 
the vector of modified residues [RP] using relations (85), (87), and (29). Solve (67) 
for [AG]. 

ggg. Update voltage angles using, [9] = [9] + [AG], and update PQ-node voltage 

magnitudes using [V] = [V] + [AV]. 
hhh* Set voltage magnitudes of PV-nodes equal to the specified values, and Increment 

the iteration count ITR=rrR+l , and proceed to step (e) 
m. Calculate line flows and output the desired results 
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Four lettered steps are characteristic steps of algorithni-4. This method is useful particularly 
for distribution systems without PV-nodes. Fig.4 is the flow-chart embodiment of 
algorifhm-4. 

Common Statements Concemb^ all methods: 

In all the prior art and invented models [¥6] and [YV] are real, sparse, symmetrical and 
built only from network elements. Since they are constant, they need to be factorized once 
only at the start of the solution. Equations (1) and (2) are to be solved repeatedly by 
forward and backward substitutions. 

[Y9] and [YV] are of the same dimensions (m+k) x (m+k) when only a row/column of the 
slack-node is excluded and both are triangularized using the same ordering regardless of 
the node-types. For a row/column corresponding to a PV-node excluded in [YV], use a 
large diagonal to mask out the effects of the off-diagonal terms. When the node is switched 
to the PQ-type the row/column is reactivated by removing the large diagonal. This 
technique is especially useful in the treatment of PV-nodes in the mattix [YV]. 

It is invented to make this technique efficient while solving (5) or (65) for [AV] by 
skipping all PV-nodes and factor elements with indices corresponding to PV-nodes. In 
other words efiEidency can be realized by skipping operations on rows/coluimis 
corresponding to PV-nodes in the forward-backward solution of (5) or (65). This has been 
implemented and time saving of about 4% of the total solution time (including 
input/output) could be realized in 14-14 iterations required to solve 1 18-node system with 
the uniform R-scale factor of 4 applied It shoidd be noted that the same indexing and 
addressing information can be nsed for the storage of both matrices as they are of die 
same dimension and sparsity structure. 

ALGORITHMS using GLOBAL CORRECTIONS 

The algorifhms-l, -2, -3, and -4 in the above involve incremental (or local) corrections. All 
the above algorithms can be oiganized to produce corrections to the initial estimate 
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solution. It iavolves storage of tiie vectcss of modified residues and replacing the relations 
(17), (18), (19) by (92), (93), (94) respectively, and (4) or (68) and (6) or (69) respectively 
by (95) and (96). Superscript '0' in relations (95) and (96) iadicates the initial solution 



estimate. 

RP/ = [(APpO' + (GppV Bpp') (AQpO'lW + ^v^'^ (92) 

RQp' = [(AQp")' - (GppV Bpp') (APpO'yW + RQp^"'^ (93) 

RPp'= APpV[Kp(VpO^] + RPp^^''^ (94) 

V=9p° + Aep' (95) 

V=Vp° + AVp' (96) 



RECTANGULAR COORDINATE FORMULATIONS OF INVENTED 
LOADFLOW METHODS 

This involves following changes in the equadons describing the loadflow models 
fonnulated in polar coordinates. 

(i) Replace 8 and AS respectively by f and Af in equations (1),(3), (4), (67), (68) and (95) 

(ii) Replace V and AV respectively by e and Ae in equations (2), (5), (6), (65), (66), (69) 
and (96) 

(iii) Replace Vp by Cp and V, by e, in equations (17) to (19), (22), (30), (31), (33), (34), 
(38) to (41), (45) to (47), (49) to (51), (54), (55), (57) to (59), (62), (63), (70), (71), 
(74) to (76), (79) to (81), (84) to (87), (90), (91). The subscript 's' indicates the slack- 
node variable. 
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(iv) After calculation of corrections to the imaginary part of complex voltage (Af) of 
PV-nodes and updating the imaginary component (f) of PV-nodes, calculate real 
component by: 



= -nIVsp' - fp' (97) 



APPENDIX 



Transformation of Branch Admittance 



The branch admittance transformation for symmetrical gain matrices of the methods 
described in the above is given by the following steps: 

1. Compute: Op = arctan (Gpp/Bpp) and 

<E>q = arctan (Gqq/Bqq) (98) 

2. Compute the average of rotations at the terminal nodes (p and q) of a branch: 

®av=(Op + Oq)/2 (99) 

3. Compare Oav with the Limiting Rotation Angle (LRA) and let Oav to be the smaller 
of the two: 



Oav = minimum (Oav, LRA) ( 1 00) 



4. Compute transformed pq-th element of the admittance matrix: 



Gpq' + jBpq' = (Cos Oav + j Siu <E>av) (C3pq + jBpq) (101) 
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5. Note that the transformed branch reactance is: 



Xpq' = BpqV(Gpq'^ + Bpq'^) and similarly, (102) 



Xpp' = BppV(Gp'+Bpp'') (103) 



In the description above X^q is the transfortned branch reactance defined by equation (1 03) 
and Bpq' is the corresponding transformed element of the susceptance matrix. Gpp' and Bpp* 
are diagonal elements obtained £nom (102). 

SOME POSSIBLE SIMPLE VARIATIONS OF SSDL-METHODS 

1. Simple obvious modifications are the use of Vp and Vp^ interchangeably in all 
expressions of RPp, and the use of 1.0 for Vs^ in all expressions of bp' involving the 
temi Vs^ 

2. bp' can also take values without transformation of bp and QSHp 

3. Explicit algorithmic steps are not given for many variants of SSDL-X*X' except 
SSDL-YY, They are obvious from tiieir descriptions and are similar to those of 
SSDL-YY method 



EXPLANATORY STATEMENTS 

The system stores a representation of the reactive capability characteristic of each machine 
and these characteristics act as constraints on the reactive power, which can be calculated 
for each machine. 



While the description above refers to particular embodiments of the present invention, it 
will be understood that many modifications may be ntiade without departing firom the spirit 
thereof The accompanying claims are intended to cover such modifications as woxxld fall 
within the true scope and spirit of the present invaition. 
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The presently disclosed embodimeats are therefore to be considered in all respect as 
illustrative and not restrictive, the scope of the invention being indicated by the appended 
claims in addition to the foregoing description, and all changes which come within the 
meaning and range of equivalency of the claims are therefore intended to be embraced 
herein. 
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The present invention is applicable to systems to process Loadflow computation by means 
of modified real and reactive power residues, and gain matrices derived from the Jacobian 
matrix. The embodiments of the invention in which an exclusive property or privilege is 
claimed are defined as follows: 

1 . A method of controlling security (over load, under/over voltage) in a power system, 
comprising the steps of: 

obtaining on-line data of nodal injections, voltages and phases at main nodes, and 
open/close status of circuit breakers in the power systen^i, 

establishing initial specifications of controlled parameters (real and reactive power 
at PQ-nodes, real power and voltage magnitude at PV-nodes, and 
transformer turns ratios etc.), 

perforaiing Loadflow computation at said nodes of tiie power system by a Super 
Super Decoupled computation in any of the Super Super Decoupled 
Loadflow methods or any of their hybrid combination or simple variants 
employing corresponding gain matrices derived from a super decoupled 
Jacobian matrix for real power with respect to angle and a sup^ decoupled 
Jacobian matrix for reactive power with respect to voltage, and involving 
triangular frtctorization of said gain matrices and computing of discrepancy 
of real power and reactive power from specified values through such nodes, 

said computing including introducing variables representing quotients of flie 
transformed discrepancies from specified values of real and reactive power 
flowing in through each node divided by voltage, or square of the voltage in 
case of transformed real power mismatches in methods employing (10, IV) 
iteration scheme, on each node, and using such variables to calculate values 
of angle and voltage for said transfomied discrepancies from specified 
values of real and reactive power flowing in through each node, by using 
triangular factorization of said gain matrices for real and reactive power, 
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initiating said Loadflow computation with guess solution of the same voltage 
magnitude and angle as those of the slack (reference) node referred to as 
slack start, 

restricting nodal transformation angle to maximum -48 degrees, applied to complex 
power iajection in conoputing said transformed discrepancies from specified 
values of real and reactive power flowing in through each node, 
evaluating the computed Loadflow for security (over load, under/over voltage), 
correcting one or more controlled parameters with said correction (amount of over 
load and/or undei/over voltage) values and repeating the computing and 
evaluating steps until evaluating step finds a good power system (no over 
load, no under/over voltage), and 
effecting a change in the voltages and phases at said nodes of the power system by 
actually implementing the finally obtained values of controlled parameters 
after evaluating step finds a good power systemi 

A method as defined in claiml wherein said Super Supor Decoupled methods, 
employing successive (19, IV) iteration scheme, of Loadflow con5)utation 
are characterized in modifying the transformed real power residue at a PQ- 
node and real power residue at a PV-node by dividing them by squared 
voltage magnitude multipHed by a factor determined as a ratio of a diagonal 
element of susceptance matrix to a diagonal element of corresponding said 
gain matrix derived from transformed Jacobian matrix for real power with 
respect to angle. 

A method as defined in claiml wherein said Super Super Decoupled methods, 
employing simultaneous (IV, 19) iteration scheme, of Loadflow 
computation are characterized in that they involve only one time calculation 
of real and reactive power residues in an iteration unlike two calculations in 
successive (19, IV) iteration scheme, and thereby achieving consequent 
speed-up. 

A simple system/method of controlling generator and transformer voltages of more 
elaborate method of security control defined in claim 1 can be realised in a 
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system for controlling generator and transfonna: voltages in an electrical 
power utility containing plurality of electromechanical rotating machines, 
transformers and electrical loads connected in a network, each machine 
having a reactive power characteristic and excitation element which is 
controllable for adjusting the reactive power generated or absorbed by the 
machine, and some of the transformers having controllable taps for adjusting 
terminal voltage, said system comprising: 

means defining any of Super Super Decoupled models of the network for providing 
an indication of the quantity of reactive power to be supplied by generators 
including at a reference node in dependence on representations of selected 
network electrical parameters, 

machine control means connected to the said excitation element of at least one of the 
rotating machines for controlling the operation of the excitation element of 
at least one machine to produce or absorb the amount of reactive power 
indicated by said model means with respect to the set of representations. 

5. A system as defined in claim 4 wherein the network includes a plurality of nodes 

each connected to at least one o£ a reference generator, a rotating machine; 
and an electrical load, and the model has one of the 3-forms and their 
variants of Super Super Decoupled matrices which receives representations 
of selected values of the real and reactive power flow firom each machine 
and to each load, and the model is operative for producing calculated values 
for the reactive power quantity to be produced or absorbed by each machine. 

6. A system as defined in claim 4 wherein the utility fiirther has at least one transformer 

having an adjustable transfom^iation ratio, and said meanis defining a model 
is further operative for producing a calculated value for the transformer 
transformation ratio. 

7. A system as defined in claim 4 wherein said machine control means are connected to 

said excitation element of each machine for controlling the operation of the 
excitation element of each machine. 
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A system as defined in claim 4 wherein said transfonnation ratio control means are 
connected to said transfomier tap changing element of each transformer for 
controlling the operation of the transfomier tap changing element of each 
transformer. 

A method for controlling generator and transformer voltages in an electrical power 
utility containing plurality of electromechanical rotating machines, 
transformers and electrical loads connected in a network, each machine 
having a reactive power characteristic and excitation element which is 
controllable for adjusting the reactive power generated or absorbed by the 
machine, and some of the transformers having controllable taps for adjusting 
terminal voltage, said method comprising: 

creating any of said Super Super Decoupled models of the network for providing an 
indication of the quantity of reactive power to be supplied by the generators 
in dependence on representations of selected network electrical parameters, 

controlling the operation of the excitation element of at least one machine to 
produce or absorb the amount of reactive power indicated by any of the said 
Super Supra: decoupled models with respect to the set of specified 
parameters. 

A method as defined in claim 9 wherein said step of controlling is carried out to 

control the excitation element of each machine. 
A method as defined in claim 9 wherein said step of controlling is carried out to 

control the tap-changing element of each transformer. 
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Fig.l : PRIOR ART: Flow-chart of FSDL solution algorithin-l 
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Fig«2: Invention: Flow-chart of SSDL-YY solution algoritiini-2 
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r START J 



Read system data and assign an initial approx. solution, or use slack-start 



Form Nodal Admittance Matrix & Initialize rnU^ITRQ^rO 



ddd 




Cos#p=Cos{-10to-20) & 
Sin$p='Sin (-10 to -20) 



Form admittance matrices \Yff] & [YV], and fectorize them as stated in step-ddd of algorithm-3 



CompfutB APi'& AQi', M*erc i=l,.., (nrt-k) and j=l,..^ 



iff 




yes 



Compute [RQJ using eqn. (70) or (75) or (80) for PQ-nodes and solve (65) 
While solving (65), slop all (he rows and colimms conesponding to PV-nodes 
Conq)ute [RP] usmg ecpi. (71) or (76) or (81) and solve (67) 



ggg 



I 



update PQ-node magnitudes using [V] = [V] + [AV] 
update PQ-node and PV-node angles using [V] = [V] + [AVj 



I 



hhh 



Set voltage magnitudes of PVntodes equal to the specified 
values, and Increment ITR»ITRP-)-l 



Calculate tine flows 
and 

output desired results 



STOP ^ 



Fig.3: Invention: Flow-chart of SSDLr-BGX' solution algoritbni-3 
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( START J 



al Read system data and assign an initial approx. solutioQ, or use slack-start 



Fonn Nodal Admittance Matrix & Initialize rrRP=rrRQ=p»K) 



Cqs^p=Cos (say -20) & 
$in$p^in(say-20) 




Foim admittance matrices {Y$] & [YV], and fiictorize them as stated in step-dddd of aIgorithm-4 

1 " 



^ Compute AP^& AQ', where f«l (nrMc) and j=l ,^,m 



Is 

AP^<e&AQj'<e7i where 
M^.,(xn+k) and j=l ,..,m 



yes 



no 



■ ■ 



Compute (RQ] using eqn. (86) for PQ^odcs and solve (65) 
While solving (65), skip all the rows and columns corresponding to PV-^iodes 
Compute [RP] using eqn. (85), (87) and (29), and solve (67) 



ggg 



update PQ-node magnitudes using [V] = [V] + [AV] 
update PQ-node and PV-node angles using [V] = [Vl + t^V] 



hhh 



Set voltage magnitudes of PV-nodes ecpial to the specified 
values, and bcrement flKprnUH-l 



m 



Calculate line flows 
and' 

output desired results 



STOP 3 



Fig.4: Invention: Flow-chart of SSDL-X'GpvX* solution aIgoritiim-4 
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turns ratio etc.) 

' 1 



30 
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Load-Flow 
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Power System State 




.50 ' 
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Evaluation, Good? 
(No over load & No 
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70 
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Correction in One or 
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Parameters 




Effect Actual Control 
To set finally obtained 
Controlled Parameters 



c 



END 
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Fig. 5: Load-Flow Computatioii in Security Control of Electrical Power System 
(PRIOR ART) 



5/6 



wo 2004/023622 



PCT/CA2003/001312 



10 

Impedance 
Network 




12 

Breaker 




14 

Admittance 




30 

Constant 


► 


Status 


— 


Matrix 


► 


Super 














Decoupled 


a 


20 


Assign TransfomGiation 


"h 


Gain Matrices, 






Ratio 







40 



TieMW 



46 



42 



Generator Loads 




62 



Machine 
Excitation Control 




Transforaier 
Tap Control 




Adjust Node 




MW 


i 


Distribution 





60 



Exercise Fast 
Super Decoupled 
Load Flow 
Algorithm 



Load Flow Solution including 
Generators Reactive Power and 
Transformers Turns Ratio 
Indications 



Fig. 6: Load-Flow Computation in Voltage Control of Electrical Power System 
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SUPER DECOUPLED LOADFLOW HETODS 

ABSTRACT 

The complete patent specification involves the invention of 
six super decoupled loadflov methods. The Fast Super Decoupled 
Loadflov (FSDL) method and two variants Transformed Fast 
Decoupled Loadflov (TFDLXB and TFDLBX) methods are known ones and 
included are their novel versions (NFSDL, NTFDLXB, NTFDLBX) • 
These are the best versions of many simple variants with almost 
similar performance. These methods specifically involved the use 
of the following invented technics in the prior art method. 

1. Gain matrices of all the six methods are different and they 
can be determined as described in the specification. 
Obviously the gain matrices of pe-subproblem can be defined 
unsymmetrical • 

2. In all the six methods, rotation angles are restricted to 
the maximum value of -36 degrees from nonlinearity 
considerations. 

3. Slack-Start for any decoupled loadflow method for efficiency 

4. Modification of real power mismatches at PV-nodes according 
to relations (12) and (15) of the specifications for the 

methods FSDL, TFDLXB, TFDLBX. 

5. - Computation of angle correction and updating for PV-nodes 

along with the voltage magnitude corrections at PQ-nodes for 
NFSDL, NTFDLXB and NTFDLBX methods. In the back-substitution 
part of the solution of this subproblem involving angle 
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corrections at PV-nodes and voltage magnitude corrections at 
PQ-nodes, skip factor elements corresponding to rows and 
columns of PV-nodes in the calculation of voltage magnitude 
corrections. 

The specification deals only with unadjusted super decoupled 
loadflow solution methods and does not describe any applications. 
However the computer algorithms can be appropriately modified 
for adjustments and/or applications such as state estimation, 
contingency analysis etc. 

Also described and claimed is the invention of two compact 
storage schemes for gain xoatrices of the FSDL method and two 
efficient procedures for node type switching implementations in 
the FSDL method. 
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THE CANADIAN PATENT ACT 



THE COMPLETE SPECIFICATION 



SUPER DECOUPLED LOADFLOW METHODS 
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The following specification particularly describes and ascertains 
the nature of this invention and the manner in which it is to be 
performed. 
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This invention relates to the super decoupled methods of 
steady-state loadflow analysis of power system, and computer 
algorithms for carrying out these methods. 

A computer algorithm processes raw information to yield 
useful information. A chemical process processes raw material to 
yield useful material. Useful information is well recognized as 
technological product in the modern age of information 
technology. Therefore, the computer algorithm like chemical 
process yields useful product and it is the useful art. 

INTRODUCTION 

Loadflow in power system studies is the most basic frequen- 
tly performed steady state analysis of an electrical power net-* 
work. The loadflows are performed in system planning, operational 
planning, and operation control. Of various methods proposed over 
the past four decades, Stott's Fast Decoupled Loadflow (FDL) [1] 
has gained wide acceptance both for off- arid on- line applica- 
tions. However, it is known to suffer poor convergence for 
systems having high ratio branches (discussions in [l] & 

[2]). 

Dy Liacco and Ramarao (discussion of [2]) as well as Deckmann 
et. al. [3] . have developed circuit transformation schemes to 
avoid such difficulties. The two schemes however do not provide 
consistent improvement in convergence [ 4 ] • Also they are not 
general and are suitable where the system has a small number of. 
troublesome R/X ratio branches. 

A recently introduced super decoupled approach [5,6,7] is 
more general and reliable. Rotation operators applied to the 
complex node injections and the corresponding admittance values 

■2 
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that relate the above to the system state variables, transform 
the network equations such that branch admittances appear to be 
almost entirely reactive. Thus, better decoupling is realized. 
However, the super decoupled algorithms are not efficient and are 
not very convenient in contingency analysis. A recent modifica* 
tion to the FDL method for networks with high R/X ratio branches 
[4] is in fact variation of the super decoupled approach with a 
weakness that it takes higher nuiober of iterations for normal 
cases. 

General-purpose version of the Fast Decoupled Loadflow 
(GFDL) method[81 proposes simple modifications to the classical 
FDL method. The convergence is much improved in the presence of 
large R/X ratio lines. The two versions of FDL method are put in 
a better theoretical framework by Monticelli et. al.[9]. The GFDL 
method is nothing but an experimental investigation of the 
original observation made about the behaviour of the FDL method 
by this inventor as back as 1985 in an unaccepted research paper 
communicated to IEEE (New York). 

The critically coupled loadflow methods are based on the use 
of the overlap update rule [10]. The disadvantage of the methods 
is that they involve the solution of a larger number of equa- 
tions at each iteration. Moreover the tests conducted for the 
paper [10] avoided the increase of R/X ratio of branches connected 
to the FV-nodes. Therefore, conclusions of the paper are not 
based on exhautive testing. Test show that the invented six 
versions are the best loadflow methods. However, simply modified 
versions including hybrids of the six invented methods can also 
have closely similar performance for any given system. 
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The latest Novel Decoupled Loadflow Method [11] can simply be 
derived from publishedtl2] and lately claimed invetion of this 
patent specification. In fact the idea involved in the method of 
[11] was revealed by this inventor in extensive communication for 
about four years with Institution of Electrical Engineers of 
England. Also the idea was investigated along with the modifica- 
tion of the real power mismatches at PV-nodes and only the best 
method was reported in reference [12]. Moreover the unsymmetrical 
gain matrix definition for the Pesubproblem is obvious from the 
definition of the factor (K) used to modify real power mismatches 
of the PV-nodes [ 12 1 . However the method of [11] is not reliable 
for systems having PV-nodes and it has no flexibility of tuning 
rotation angle as was originally known to this inventor. Possibly 
this may be the reason that the author of [11] did not comapre 
his method with the Fast Super Decoupled Loadflow [12] and 
Transformation based Fast Decoupled Loadflow [13] methods. 

Keywords : Power systems, Loadflow, Matrices, solution of 

simultaneous equations. 

The invention will now be described, wherein the following 
symbols are used ; 

y ssQ +jB : (p-q)-th element of nodal admittance 
P<3 P9 P<I matrix formed excluding shunts 

y B g +jb : total shunt admittance at node p 
P P P 

V -e +jf »V ; complex voltage at node p 
P P P P P 

IS!^ 9 ty * voltage angle, magnitude corrections 
P P 
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p p 



real, iineginary component corrections of voltage 



P + jQ 
P P 

tP +j hQ 
P P 

RP +jRQ 
P P 



: net nodal injected power calculated 



: nodal power residue (misinatch) 



: modified nodal power residue 



in 
k 

n=ro+k+l 
[ 1 



nu 



(er of PQ-nodes 



: number of PV-nodes 

: total number of system nodes 

: q is the node adjacent to node p excluding the case q 

: indicates enclosed variables to be vector or matrix 



SUPER DECpUPIiED LOADFLOW (SDL) METHOD { The prior art) 

There are two versions of the Super Decoupled Loadflow 
method. These are based on the XB- and BX- versions of the Fast 
Decoupled Loadflow (FDL) method. The one based on the XB-version 
was developed in 1985 [7] and the other based on the BX-version 
is obvious from the BX-version of the FDL method developed in 
1989 [81. These methods involve the iterative solution of system 
of equations (1) and (2). Unlike FDL model the gain matrices in 
its. transformed version ^ are unsymmetrical because mostly 
different rotations are required to be applied at the terminal 
nodes of a branch. Symmetrical gain matrices can be obtained by 
the Haley and Ayres technique [7] of applying average of the 
rotations at the terminal nodes of a branch to the branch 
admittance (appendix) • 
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SDL based on the XB-version of the FDL (SDLXB) 

The SDLXB*verslon is described by the following equations 
(1) to (10). 

[RP] = [Y'] [ t^e] (1) 

[RQ] = [Y"] [ (2) 

Wherein each element of [RP] and [RQ] at PQ-nodes are given by 
(3) and (4) respectively. Whereas an element of [RP] at FV-nodes 
is given by (5) • 

RP = ( fciP Cosp + ^^Q Sinp )/V (3) 
P P P P P P 

RQ = ("h^ Sinp + ^Q Cosfi )/V (4) 
P P P P P P 

RP = ts.P /V (5) 
P P P 

Trigonometric functions in (3) and (4) and elements of [Y'] and 



[Y"] are given by (6), (7) and 


, (8), (9) and (10). 




f 2 2 
COBp = -B / /g + B 

P PP V PP PP 




(6) 


r2 2 
Sinfi « -G / /G + B 
P PP V PP PP 




(7) 


Y' = -1 / X and 

pq pq 


Y" = -B 

pq pq 


(8) 


Y' =Z"^' 

PP T-^f pq 


Y" = -2b '+ X -Y" 

pq p *f*r pq 




b' = b Cos^ or = b 
P P P P 


• 


(10) 


Where X is the transformed branch reactance defined in 

pq ^ 

by eqation (41) and B is the corresponding transformed 


appendix 
element 



pq. 

of the susceptance matrix. 



6 
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SDL based on the BX-version of the FDL (SDLBX) 



The SDLBX-version differs from the SDLXB-version only in 
relation (8) as given by relation (11). The SDLBX-version 
consists of relations (1) to (7), (11), (9) and (10) • 



-B 



and 



pq 



pq 



Y" = -1 / X 

pq pq 



(11) 



INVENTED SUPER DECOUPLED LOADFLOW METHODS 

Six invented super decoupled loadflow models are described 
in this section. Fast Super Decoupled and two versions of *the 
Transformed Fast Decoupled models are Known ones [12,13] and 
included are their novel versions. 
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Fast Super Decuopled Loadflow (FSDL) 



The FSDL method involves the iterative solution of the 
system of equations (1) and (2). The model cosists of relations 
(1) to (4), (6), (7), (12) to (16), (9) and (10). 



*^ ■ • » V 



• :• •■ .-r 



RP 



/(K V ) 
P P P 



(12) 



Elements of [Y'J and [Y"]. and the multiplier K in (12) are given 

P 

by : 



pq 



-Y 



pq 



for branch r/x ratio < 2.0 



-{B +0.9(Y -B )) for branch r/x ratio > 2.0 

pq pq pq 



-3pq B 



for branches connected between 
two PV-nodes or a PV-node and 



the slack-node 



(13) 



- v.* 



Y" 



pq 




-Y 



pq 



for branch r/x ratio $ 2.0 



-(B +0.9(Y -B )) for branch r/x ratio > 2.0 
pq pq pq 

7 



(14) 
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K = Absolute (B /Y' ) (15) 
P . PP PP 

KKp = Absolute (Y'pp/Y"pp) (16) 



Branch admittance magnitude in (13) and (14) is of the sane 
algebraic sign as its susceptance. Elements of the two gain 
matrices differ in that diagonal elements of [Y"] additionally 
contain the b' values given by equation (10) and in respect of 
elements corresponding to brcuiches connected between two PV*nodes 
or a FV-node and the slack-^node. The factor 0.9 in the relations 
(13) and (14) is tuned only for rotation limited to -36 degrees. 
With different Limiting Rotation Angle (LRA), it needs to be tuned 
again. In two simple variations of the FSDL method, one is to 
make Y" = Y' and the other is to make Y' » Y" 

pq pq pq pq 

Transformation based Fast Decoupled Loadf low (TFDL) 

The TFDL model is similar to the FSDL model. They differ 
only in the definition of the gain matrices. 
The TFDL(XB-version) (TFDLXB) Loadflow 

The TFDLXB-version consists of relations (1) to (4), (12), 
(15), (6), (7), (17), (18), (19), (9) and (10). 



-Bpg for branches connected between two FV- 

nodes or a PV-node and the slack-node 



Y'pq = 

-1 / X for all other branches (17) 

— pq 

Y" « -B (18) 

pq pq 

KK = Absolute (Y' /(51l/X )) (19) 

p PP pq 

8 
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The TrDL(BX-version) (TFDLBX) Loadflow 



The TFDLBX-version consists of relations (l) to (4), (12), 
(15), (6), (7), (20), (21), (22), (9) and (10). 



. ■ * ■ • - • » ■ - 



Y'pq = 



-Bpq 



B 



pq 



pq 



-1 / X 



pq 



for branches connected between two PV- 
nodes or a PV-node and the slack-node 



for all other branches 



(20) 



(21) 



* • • .» ^» • * • , 
^ " • ' ••*. * • 

• . • ■ ... I..' 1 



KK = Absolute (Y' /(21b )) 
P PP 1-»f pq 



(22) 



Where X is the transformed branch reactance defined in appendix 

pq 

by eqation (41) and B is the corresponding transformed element 

pq 

of the susceptance matrix* 
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The factor KK is to be multiplied to the real power mis- 

P 

match at a PV-node switched to PQ-type in node type switching 

implementation. Prom general considerations, K and KK are 

P P 

restricted to the minimum value of 0.75 for FSDL, TFDLXB and 
TFDLBX methods. However they can be tuned for the best possible 
convergence for any given system. 



. ; « ■•J*, * B ■ • 



. • / . • * • 



' •'. .'- •-. 



! I 



In all the SDLXB, SDLBX, FSDL, TFDLXB and TFDLBX models 
[Y'J and [Y"] are real, sparse, symmetrical and built only from 
network elements. Since they are constant, they need to be 
factorized once only at the start of the solution. Equations (1) 
and (2) are to be solved repeatedly by forward and backward 
substitutions • 
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[Y'] and [V"] are of the same dimensions (m-fk}*(m+k) when 
only a row/column of the slack-node is excluded and both are 
triangularized using the same ordering regardless of the node- 
types. For a row/ column corresponding to a PV-node excluded in 
[Y"], use a large diagonal to mask out the effects of the off- 
diagonal texnns. When the node is switched to the PQ-state the 
row/column is reactivated by removing the large diagonal* This 
technique is especially useful in the treatment of PV-nodes in 
the matrix [Y"]. 

It is invented to make this technique efficient while 
solving (2) for ^V by skipping all PV-nodes and factor elements 
with indices corresponding to PV-nodes. In other words efficiency 
can be realized by skipping operations on rows/columns 
corresponding to PV-nodes in the forward-backward solution of (2) 
for ^V. This has been implemented and the time saving of about 
4% of the total solution time (including input/ output) could be 
realized in 14-14 iterations required to solve 118-node system 
with the uniform R-scale factor 4 applied. The time saving has 
been assessed on PC-XT. It should be noted that the same indexing 
and addressing information can be used for the storage of both 
the matrices as they are of the same dimensions and sparsity 
structure • 



Novel Fast Super Decoupled Loadflow (NFSDL) 

The NFSDL method involves the iterative solution of the 
following system of equations (23) and (24). 

[RPJ = [Yl] [ tfii (23) 



ERQ(PQ-nodes7l = rY2l fty^V ( PQ-nodes )1 (24) 
RP(PV-nodes2J L J |^©(PV-nodes2j 
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Wherein each element of [RP] and [RQ] at PQ-nodes are given by 
(3) and (4) respectively. Whereas an element of [RP] at PV^nodes 
is given by (5)* In (24) PV-nodes are assumed to be grouped and 
numbered after all PQ-nodes. 
Elements of [Yl] and [Y21 are given by : 



pq 



for branch r/x ratio ^2.0 



pq 



-(B +0.9 (Y -B )} for branch r/x ratio > 2.0 



pq 



-Bpq 



pq pq 



for branches connected between 
two PV-nodes or a PV-^node and 
the slack-node (25) 



Y2pq = 



Y2pp 



Gpq 




Yl 



for branches connected to a PQ-node 
and a PV-node 



-Ylpq for branches connecting two PQ-nodes 



for branches connecting two PV-nodes 



pq 



Yl =21 -Yl 
pp <i^p pq 



2bj^-Ylpp 



for PQ-nodes 




for PV-nodes 



(26) 



(27) 



(28) 



-B 

- PP 

The NFSDL model consists of relations (23), (24), (3) to 
(7), (25) to (28), and (10). Branch admittance magnitude in 
(25) is of the s£une algebraic sign as its susceptance. 



Novel Transformed Fast Decoupled Loadf low (NTFDL) 



The NTFDL model is similar to the NFSDL model. They differ 

only in the defination of gain matrices. 

11 
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The NTFDL(XB-version) (NTFDLXB) Loadflow 

The NTFDLXB model consists of relations (23), (24), (3) to 
(7), (29) to (32), and (10). 



Ylpq = 



-Bpq 



-1 / X 
L- pq 



for branches connected between two PV- 
nodes or a PV-node and the slack-node 



for all other branches 



(29) 



Y2pq =[Gpq 



Yl 



PP 



B 



pq 



-B 



= S-Y1 

pq 



for branches connected to a PQ-node 
and a PV-node 



for branches connecting two PQ-nodes 



for branches connecting two PV-nodes (30) 



(31) 



Y2 



PP 



2b '+B 
P PP 



-B 



PP 



for PQ-nodes 



for PV-nodes 



(32) 



Where X Is the transformed branch reactance defined in appendix 

pq 



by eqation (41) and B is the corresponding transformed element 

pq 

of the susceptance matrix. 

The NTFDL(BX-version) (NTFDLBX) Loadflow 

The NTFDLBX model consists of relations (23), (24)^ (3) to 
(7), (33) to (36), and (10). 




%•• V....;! 
'.•..-.•/.••.viy, 



Ylpq = 



-Bpq 



for branches connected between two PV- 
nodes or a PV-node and the slack-node 



.»'.•.-••-.••,■••'.•.• 



-B 



for all other branches 



L 



pq 



(33) 
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Y2pq 



Yl 



PP 



Y2 



PP 



Gpq 



1/ X 



pq 



-B 

Lpq 



pq 



for branches connected to a PQ-node 
and a PV-node 

for branches connecting two PQ-nodes 



for branches connecting two PV-nodes (34) 



2b' +Z.-V X 

p t->i> pq 



or « 2b'+l/X for PQ-nodes 
P PP 



(35) 



-B 



PP 



for PV-nodes (36) 



In the NFSDL NTFDLXB and NTFDLBX methods [Yl] and [Y2] are 
real, sparse, symmetrical and built only from network elements. 
Since they are constant, they need to be f actor ized once only at 
the start of the solution. Equations (23) and (24) are to be 
solved repeatedly by forward and backward substitutions. 

[Yl] and [Y2] are of the same dimensions (m+k) * (m+k) when 
only a row/column of the slack-node is excluded and both are 
triangularized using the same ordering regardless of the node* 
types. It should be noted that the same indexing and addressing 
information can be used for the storage of both the matrices as 
they are of the same dimensions and spars ity structure. Unlike 
[Y"}, It is to be noted that all the PV-nodes are also active in 
[Y2]. Therefore, the novel methods would pose a difficult problem 
of .node-type switching. 

p is restricted to the maximum of -36 degrees from 
P 

nonlinear ity considerations for FSDL, TFDLXB, TFDLBX, NFSDI*, 
NTFDLXB and NTFDLBX methods. However it can be tuned for the 
best possible convergence for any given system and it can also be 



vertually unrstricted -90 degrees. 

. . 13 
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Iteration Scheme 

The basic iterative scheme for solving the proposed models 
is to .solve for [^6] to update [6] and then solve (2) or (24). 
This is the block Gauss-Seidel approach. The scheme is block- 
successive which imparts increased stability to the solution 
process. This in turn improves convergence and increases 
reliability. Cycling behaviour in the iterative process of the 
standared iteration scheme of [1] can be avoided by simply using 
different convergence tolerances for the real and reactive power 
mismatches as suggested in [9]. The cycling behaviour and its 
remedy as above were originally observed by this inventor as back 
as in 1985 in an unaccepted paper submitted to IEEE and prepared 
from research conducted at the university of Roorkee* However, 
algorithms are given for strictly successive iteration scheme of 
[8]. 

APPENDIX 



2. Compute the average of rotations at the terminal nodes (p 
and g)of a branch : 



* • * - 



Transformation of Branch Admittance 

The branch admittance trans foramat ion for symmetrical gain 
matrices of the TFDLXB, TFDLBX, NTFDLXB and NTFDLBX methods is 
given by the following steps : -^i-^S 



.". y. . ..V..-.V 



1. Compute : 0 « arctan (G /B ) and 

PP PP PP 

(37) 

p = arctan (G /B ) 

qq qq qq -vigt 



"V. . 



P = + P ) / 2 (38) 

av PP qq 
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3. Compare p with the Limiting Rotation Angle (LRA) and let 

av 

to be the smaller of the two : 

av 

fi - minimum [fi , LRA) (39) 
av av 

4. Compute transformed pq-th element of the admittance matrix : 

6 +jB = (Cosp +jSinJJ ) (G +jB ) (40) 
pq pq av av pq pq 

5. Note that the transformed branch reactance is : 

2 2 

X = B / (G +B ) (41) 
pq pq pq pq 

and similarly 

2 2 

X = B / (G +B ) (42) 
PP PP PP PP 



ALGORITHM-1 (The prior art) 

Super Decoupled Loadflow Algorithms (solution Steps) 
Strictly Successive (18, IV) Iterative Scheme : 

(a) Read system data and assign an initial approximate solution. 
If better solution estimate is not available, set specified 
voltage magnitudes at PV^nodes and 1.0 p.u. voltage 
magnitudes at PQ-nodes. Set all the node angles equal to 
that of the Slack-node angle. This is referred to as the 
flat-start [1]. 

(b) Initialize iteration counts ITRP=ITRQ=r=0 

(c) Form nodal admittance matrix 

15 
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(d) Form (m+k) by (m+k) size matrices [Y'J and [Y"] of (1) and 
(2) respectively each in a compact storage exploiting 
sparsity. 

(i) In case of SDLXB-method, the matrices are formed usig 
relations (8), (9) and (10) • 

(ii) In case of SDLBX-method, the matrices are formed using 
relations (11), (9) and (10). 

(e) In [Y"] matrix, replace diagonal elements corresponding to FV 
-nodes by very large value (say 10 •0**10). In case [Y"] is of •t!!^^ 
dimension (m by m),this is not reqired to be performed. 
Factorize [Y'j and [Y"] using the same ordering regardless 
of the node-types. In case [Y**] is of dimension (m by m) , it 
is factorized using different ordering than that of [Y^]. 



nodes only) • If all are less than the tolerance ( ) / 



(g) Compute the vector of modified residues [RP] using (3) for 
PQ-nodes and using (5) for PV-nodes. 

(h) Solve (1) for voltage angle corrections and update voltage 
angles using eqn. (43) • 

« 

(i) Increament the iteration count ITRP»ITRP+1 and 
r=»(ITRP+ITRQ)/2. 

r r 
(j) Compute residues (PQ- and PV~nodes) and (at PQ- 

nodes only) . If all are less than the tolerance ( ) # 

proceed to step (n) • Otherwise follow the next step. 

16 
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(f) Compute residues ^P (PQ- and PV-nodes) and (at PQ- 



proceed to step (n) . Otherwise follow the next step. 
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(k) Compute the vector of modified residues [RQ] using (4) for 
PQ-nodes only. 

(1) Solve (2) for voltage magnitude corrections and update 
voltage magnitudes of PQ*nodes only using eqn.(44). 

(m) Increament the iteration count ITRQ=ITRQ+1 and 
r=(ITRP+ITRQ)/2. Proceed to step (f) 

(n) Calculate line flows and output the desired results. 

The SDLXB and SDLBX algorithms differ only in step-d defining 
gain matrices in the above solution steps. Fig.l is the flowchart 
of algorithm-1. 

ALGORITHM-2 (containing invented steps) 

Fast Super Decoupled Loadflow Algorithms (solution Steps) 
Strictly Successive (19, IV) Iterative Scheme | 

(a) Read system data and assign an initial approximate solution. 
If better solution estimate is not available, set all the 
nodes voltage magnitudes and angles equal to those of the 
Slack-node. This is referred to as the slack-start. 

(b) Initialize iteration counts ITRP^ITRQ-r^O 
. (c) Form nodal admittance matrix 

(d) Form (m+k) by (m+k) size matrices [Y^J and [Y"] of (1) and 
(2) respectively each in a compact storage exploiting 
spars ity. 

(i) In case of FSDL-method, the matrices are formed using 
relations (13), (14), (9) and (10). 
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(ii) In case of TFDLXB-method, the matrices are formed usig 
relations (17), (18), (9) and (10). 

(iii) In case of TFDLBX-method, the matrices are formed 
using relations (20), (21), (9) and (10). 

(e) In [Y"] matrix, replace diagonal elements corresponding to PV 
-nodes by very large value (say 10. 0**10). In case [Y"] is of 
dimension (m by m) , this is not reqired to be done. 
Factorize [Y'] and [Y"] using the same ordering regardless 
of the node-types. In case [Y»] is of dimension (m by m) , it 

•■•••■•x-'.'V:/': 

is f actor ized using different ordering than that of [Y'J* S^^S^S 



(g) Compute the vector of modified residues [RP] using (3) for 



(h) Solve (1) for voltage angle corrections and update voltage 
angles using relation (43) • 

r (r-l) r 

e = e + (iB (43) 
P P P 

(i) Set the voltage magnitudes of PV-nodes equal to the 
specified values. This is required to be done only in the 
first iteration when slack-start is used. Increament the 
iteration count ITRP=ITRP+1 and r=C^'^'*-*- J^TR^j/z- 

"i 



J -: - • 

i' « 



• •mm • • 



(f) Compute residues (PQ- and PV-nodes) and (at PQ- 

nodes only) . If all are less than the tolerance ( £^ ) , 
proceed to step (n) . Otherwise follow the next step. 




' :'.■ :<•■-•■• 
• ,•• • /••. mV 



PQ-nodes and using (12) and (15) for PV-nodes. In (3) p is 

P :^>S&> 

restricted to the maximum of -36 degrees (or any other 
angle) from nonlinearity considerations. It can be even 
unrestricted for any given system (-90 degree) . 
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r r 

(j) Compute residues (^P (PQ- and PV-nodes) and fsJQ (at PQ- 

nodes only) • If all are less than the tolerance ( £ ) / 
proceed to step (n) • Otherwise follow the next step. 

(k) Compute the vector of modified residues [RQ] using (4). In 

(4) 0 is restricted to the maximum of -36 degrees (or ^ any 
P 

other angle) from nonlinearity considerations. It can be 
even unrestricted for any given system (-90 degree) . 

(1) Solve (2) for voltage magnitude corrections and update 
voltage magnitudes of PQ-nodes only as given by relation 
(44). While solving (2), skip all the rows and columns 
corresponding to PV-nodes to increase the efficiency. 

V « V + (iV (44) 

P P P 

(m) Increament the iteration count ITRQ=ITRQ+1 and 
r=(ITRP+ITRQ)/2. Proceed to step (f) 

(n) Calculate line flows and output the desired results. 

All the above steps can also be the parts of the classical 
(le, IV) iteration scheme of reference [1]. This algorithm is for 
solving loadflow problem formulated in polar coordinates. The 
solution of the loadflow model formulated in rectangular 
coordinates also Involve similar steps with minor modifications. 
Gain matrices of the invented methods are the same for both 
polar- and rectangular-coordinate formulations of the loadflow 
problem. Fig. 2 is the flowchart of algorithm-2. The FSDL-method, 
TFDLXB-method and TFDLBX-method differ only in step-d defining 
gain matrices in the above algorithm-2. 
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ALGORITHM-3 (invclude invented steps) 

Novel Fast Super Decoupled Loadflow Algorithms (solution Steps) 
Strictly Successive (18, IV) Iterative Scheme | 

(a) Assign an initial approximate solution. If better solution 
estimate is not available, set all the nodes voltage 
magnitudes and angles equal to those of the Slack^node. This 
is referred to as the slack-start* 

(e) Initialize iteration counts lTRP=ITRQ=r=0 

(c) Form nodal admittance matrix 



(m+k) by (m+k) size matrices [Yl] and [Y2] of (23) and 
respectively each in a compact storage exploiting 



(d) Form 
(24) 
sparsity. 

(iv) In case of NFSDL-method, the matrices are formed using 
relations (25) to (28), and (10). 

(v) In case of NTFDLXB-method, the matrices are formed 
using relations (29) to (32) and (10) . 

(vi) In case of NTFDLBX-method, the matrices are formed 
using relations (33) to (36) and (10). 

(e) Factorize [Yl] and [Y2] using the same ordering regardless 
' of the node-types. 



V • 
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(f) Compute residues (^P (PQ- and PV-nodes) and (at PQ- 

nodes only). If all are less than the tolerance ( £^), 
proceed to step (n) • Otherwise follow the next step. 
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(g) Compute the vector of loodified residues [RP] using (3) for 

PQ-nodes and using (5) for PV-nodes. In (3) p is restricted 

P 

to the maximum of -36 degrees (or any other angle) from 
nonlinear ity considerations. It can be even unrestricted for 
any given system (-90 degrees) . 

(h) Solve (23) for voltage angle corrections and update voltage 
angles using relation (43). 

(i) Set the voltage magnitudes of PV-nodes equal to the 
specified values. This is required to be done only in the 
first iteration when slack-start is used. Increament the 
iteration count ITRPs=ITRP+l and r=(ITRP+ITRQ)/2. 

r r 
(j) Compute residues ^P (PQ- and PV-nodes) and (at PQ- 

nodes only) . If all are less than the tolerance ( ) r 

proceed to step (n) • Otherwise follow the next step. 

(k) Compute the vector of modified residues [RQ] using (4) for 

PQ-nodes. Compute [RP] for PV-nodes using (5).In(4)p is 

P 

restricted to the maximum value of -36 degrees (or any other 
angle) from nonlinearity considerations. It can be even 
unrestricted for any given system (-90 degrees). 

(1) Solve (24) for voltage magnitude corrections at PQ-nodes 
and voltage angle corrections at PV-nodes, Update voltage 
magnitudes of PQ-nodes using relation (44) and angles of PV- 
nodes using relation (43) . In the back substitution part ot 
the solution of (24), skip all the PV-nodes in the 
calculation of voltage-magnitude corrections of PQ-nodes. 

21 
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(m) Increament the iteration count ITRQ=ITRQ+1 and 
r=(lTRP+ITRQ)/2. Proceed to step (f) 

(n) Calculate line flows and output the desired results. 

All the above steps can also be the parts of the classical^ 
(le, IV) iteration scheme of reference [1]. This algorithm is for 
solving loadflow problem formulated in polar coordinates. The 
solution of the loadflow model formulated in rectangular 
coordinates also involve similar steps with minor modifications. 
Gain matrices of the invented methods are the same for both 
polar- and rectangular-coordinate formulations of the loadflow 
problem. The NFSDL-method, NTFDLXB-method and NTFDLBX-method 
differ only in step-d defining gain matrices in the above 
algorithm-3* Fig. 3 is the flow-chart of algorithm-3. 

ALGORITHMS using GLOBAL CORRECTIONS 

The algorithms-1, -2 and -3 of above involve increament al (or 

local) corrections. All the above algorithms can be organised to 

produce corrections to the initial estimate solution. It involves 

storage of the vectors of modified residues and replacing the 

relations (3) , (4) , (5) , (12) , (43) and (44) respectivly by 

(45) , (46) , (47) , (48), (49) and (50). Superscript '0' in relations 

(49) and (50) indicates the initial solution estimate. 

r r r r (r-1) 

RP = ( £iP Cosfi + Sinfi )/V + RP (45) 
P P P P P P P 

r r r r (r-1) 

RQ = (-^P Sin^l + fc^Q Cosfi )/V + RQ (46) 
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r r r (r-1) 

RP « KP /V + RP (47) 

P P P P 

r r r (r-l) 

RP = iiP /(K *V ) + RP (48) 

P P P P P 

r Or 

e = e + i^e (49) 

p P p 

r 0 r 

V « V + J^v (50) 

P P P 

RECTANGULAR COORDINATE FORMULATIONS OF THE INVENTED LOADFXiOW 
METHODS 



This involves following changes in the equations describing 
the loadflov models formulated in polar coordinates • 

(i) Replace e and 1^0 respectively by f and j^f in equations 

(1) , (23), (24), (43) and (49). 

(ii) Replace V and ^V respectively by e and {^e in equations 

(2) , (24), (44) and (50). 

(lii) Replace V by e or e in equations (3), (4), (5), {12), 

p p s 

(45), (46), (47) and (48). The subscript 's' indicates the 
Slack-node variable. 

(iv) After calculation of corrections to the imeginary. part of 
complex voltage ( ^f ) of PV-nodes and updating the 
imaginary component (f) of PV-nodes, calculate real 
component by : 

2 2 

e = V - f (51) 

p p( specified) p 
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COMPACT STORAGE SCHEMES FOR GAIN MATRICES OF THE FSDL METHOD 

Two compact matrix storage schemes for gain matrices of the 
FSDL method are described in the following. 

(1) When network shunts are ignored (FK=0.0), [Y"] becomes the 
submatrix of [V]. in this scheme, [Y'] is factorized using 
optimal ordering regardless of node-types along with ' the 
following procedure : 

1. Identify PQ«-nodes in the composite path of all PV-nodes« 

2. Reload and ref actor all row/columns of these PQ-nodes using 
partial refactorization method PRl of [13]. Store these 
factors separately. 

3. When solving for jj^V using factors of [Y'] : 

(i) all PV-nodes and factor elements with indices 
corresponding to PV^nodes are skipped. Zeroing of reactive 
power mismatches at PV-nodes is not required. 

(ii) for PQ-nodes of the path, separately stored factors as 
in 8t6p-2 above are used along with the other factors of 
[Y'] in forward-backward operations. 

The above procedure is implemented. For X18Tnode system, 
there are only 18 PQ-nodes in the composite path of all 53 PV- 
nodes. Each of the gain matrices [Y'] and [Y"] is having 117- 
row/columns for 118-node system, and the total of 2 34 -row/columns 
are to be factorized and stored for conventional 2 -matrix 
solution. Threfore the proposed second scheme can be said to^ be 
saving 42% (99/234='0.42) of the factorization effort and core 
memory. This saving rises to the maximum 50% for systems without 
PV-nodes . 
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(2) When network shunts are ignored (FK=»0.0}, [Y**] becomes the 
siibmatrix of [Y']. In the single matrix version of the FSDL 
method only [Y'] is factorized and stored. This is achieved 
exactly by placing PV-nodes after all PQ-nodes and by allowing 
optimal ordering only within each group. While solving (2) for 
^V, all FV-nodes and factor elements with indices corresponding 
to PV-nodes are skipped. Zeroing of reactive power mismatch at 
PV-nodes is not required. Apperantly^ 50% reduction in the matrix 
storage is achieved in the absence of PV-nodes. 

PV-node Q-limit Adjustment for the FSDL method 

When the power flow solution is moderately converged, any Q- 
limit violations should be corrected. Any violated PV-node is 
converted to the PQ-type with the rective generation set at the 
limiting value. The voltage of the converted node is subsequently 
compared to the scheduled value and the node is reconverted to 
the PV-type, if any of the back-off conditions are satisfied. 

The conventional node-type switching approach is more 
suitable for the handling of the PV-node Q-limiting. The 
switching of node-types involves the insertion/delation of 
eqalities into/from (2). This can be achieved efficiently by 
partial matrix refactorization and factor » updating methods. 
Partial Refactorization Hethod-1 (PR1)[13] can be used to update 
[^"1 in order to account for changes in the node status (regu- 
lated or nonregulated) • The efficiency of PRl is inversely 
proportional to the number of node status changes. Two efficient 
techniques of the node-type switching are described in the 
following. 
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(1) When network shunts are ignored (FK-0.0) or accounted by 
alternative ways, a reactivated row/ column by the removal of the 
large diagonal corresponding to a PV-node in [Y"] is identical 
to that in [Y'], It is this feature of the FSDL method exploited 
to enhance the efficiency of the node type switching using the 
following procedure. This is convenient with the simple variant 
of the FSDL method mentioned before* 

1. find the composite path of PV-nodes switched to PQ*type 

2. find the composite path of the rest of the PV-nodes 

3. ref actor only rows/columns of [Y"] that constitute common 
path of steps 1 and 2* The factored rows/columns of [Y'] in 
the uncommon path of step 1 can be used for the partial 
refactorization of [Y"] 

4. while solving (2), use factors from [Y'] in the uncommon 
part of the path of step 1 along with other factos of [Y"]. 

When all the PV-nodes of the system are switched to PQ-type, 
[Y'] is used for solving (2) and no partial refactorization in 
[Y"] is required. Also when common path to be determined in step- 
3 does not exist, partial refactorization is not required. 

(2) The compact storage scheme (1) for the FSDL method described 
above inherently provides for efficient implimentation of node- 
type switching used for FV-node Q-limit adjustment. The procedure 
constitutes the following steps. 

1. PQ-nodes in the composite path of all PV-nodes are known and 

their rows/columns are factored and stored separately from 

those of [Y'J for unadjusted solution. 
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2. Identify PQ-nodes in the composite path of current PV-nodes. 

3. Identify PQ-nodes which are not conuaon in steps 1 and 2* 

4. Carry out Reverse Correction Partial Matrix Refactorization 
along the composite path of PQ-nodes of step-3 on separately 
stored factors of step-1. Use large diagonals in 
rows/ columns of PQ-nodes of 8tep-3. 

If no PQ-node is identified in step-3, Reverse correction Partial 
matrix refactorization of step-4 is not required to be performed. 
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I CLAIM : 

(a) TO have invented three Super Decoupled Loadflow methods 
(FSDL, TFDLXB and TFDLBX) both in Polar coordinate and 
Rectangular coordinates formulations, their Novel versions 
(NFSDL, NTFDLXB and NTFDLBX) . and their all possible hybrid 
cobinations. Also I claim to have invented the technique of 
producing Global corrections as given by relations (45) to 
(50) . 

(b) specifically the invention involves the use of the following 
seven (7) items in the methods of (a): 

(i) the definitions of the gain matrices of the six-methods 
of (a) corresponding to steps-d in algorithm-2 and 
algorithm-3. for both polar and rectangular formulations of 
the loadflow problem. Though explicit definitions, of 
unsymmetrical gain matrices of Pe-subproblem are not given, 
they are obvious for all the six methods and they are also 
claimed by this inventor as simple modifications. 

(ii) restriction of rotation angle fi to the maximum of -36 

degrees (any other angle) in relations (3), (4), (45) and (46) 
from nonlinearity considerations, for both polar and 
rectangular coordinate formulations of the loadflow problem. 

(iii) modification of real power mismatch at PV-nodes accor- 
ding to relations (12) and (15) in FSDL, TFDLXB and TFDLBX- 
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(vii) The six invented methods of unadjusted loadflow 
solution can be used for adjusted solution by adding some 
more steps in the algorithms; state estimation, contingency 
analysis and in variety of advanced power network analysis 
and control. 
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methods formulated in both polar and rctangulr coordinates. 
This in general is the technique to modify the known vector 
to ensure the coefficient/gain matrix symmetrical. It can 
be used in all possible computer algorithms (including all 
other approaches of the formulation of the loadflow problem) 
and real-time/ on-line/off-line applications. 

(iv) computation of angle corrections and updating for PV- 
nodes along with the Voltage magnitude corrections while 
solving equation (24) for NFSDL, NTFDLXB and NTFDLBX 
methods. In the back-substitution part of the solution of 
(24), I claim to have invented the technique to skip all 
the rows and columns corresponding to PV-nodes for calcula- 
tion of voltage magnitude corrections of PQ-nodes. This also 
applies when the novel methods are formulated in rectangular 
coordinates* 

(V) the Slack-start procedure for all the decoupled loadflow 
methods and in particular for the six methods of (a) aibove. 



.'..1 ''"v >A : 



(vi) similar decoupled methods as of (a) dan also be used :^^yj^ff:- 
for solving siinultaneous equations appearing in other 
of analysis, operation and control. 
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(c) Also I claim to have Invented two compact storage schemes 
for gain matrices of the FSDL method and two efficient 
procedures of the node*type switching implementations for 
the FSDL method as described in the above. 



Dated this 7 th day of november 1993. 



Signature of the Inventor S.B.Patel 
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Fig. 1 Flow-chart of unadjusted load flow solutk>n algorlthm-1 . 
It is the same for SDLXB and SDLBX methods. (flow;:chart of the prior art) 
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Fig. 2 Flow-chart of unadjusted load flow solutton algorlthm-2. ^ 
It Is the same for FSDI-, TFDLXB and TFOIBX methods, (flowchart of Invited methods) y^c^ 
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